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ABSTRACT 

We present a high spectral resolution survey of the most metal-poor damped Lya 
absorption systems (DLAs) aimed at probing the nature and nucleosynthesis of the 
earliest generations of stars. Our survey comprises 22 systems with iron abundance 
less than 1/100 solar; observations of seven of these are reported here for the first time. 

■ Together with recent measures of the abundances of C and O in Galactic metal-poor 

stars, we reinvestigate the trend of C/O in the very metal-poor regime and we com- 

q \ pare, for the first time, the O/Fe ratios in the most metal-poor DLAs and in halo stars. 

$h ' We confirm the near-solar values of C/O in DLAs at the lowest metallicities probed, 

^ ' and find that their distribution is in agreement with that seen in Galactic halo stars. 

d . We find that the O/Fe ratio in very metal-poor (VMP) DLAs is essentially constant, 

and shows very little dispersion, with a mean [(O/Fe)] = +0.39 ±0.12 , in good agree- 
ment with the values measured in Galactic halo stars when the oxygen abundance 
is measured from the [O i] A6300 line. We speculate that such good agreement in the 
ly-j ' observed abundance trends points to a universal origin for these metals. In view of this 

\ agreement, we construct the abundance pattern for a typical very metal-poor DLA and 

OO ■ compare it to model calculations of Population II and Population III nucleosynthesis 

' to determine the origin of the metals in VMP DLAs. Our results suggest that the 

, most metal-poor DLAs may have been enriched by a generation of metal-free stars; 

■ however, given that abundance measurements are currently available for only a few 

elements, we cannot yet rule out an additional contribution from Population II stars. 
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5» , 1 INTRODUCTION erations of stars, depend on various (largely unknown) 

factors including the form of the primordial stellar ini- 
The initial conditions tor cosmic chemical evolution are of , . , r ,. , . r P . , 

tial mass function and the uniformity ol the enrichment 

fundamental importance to our understanding of galaxy ,,, . . , ,. ,. nr^Air T~- } llnnn J 

, f P- °f the mtcrgalactic medium (1GJV1; Bromm & Larson 2004; 

formation and the process of galactic chemical evolution. nr; — ; „ — — 1T~FTi — — Inn-, iK T 1 

_, ,. . , . _ , . Karlsson, Bromm, & Bland-Hawthorn 2011). In order to pm 

these conditions, set by the yields of the first few gen- , . ... , .... „ , . , . ,. 

down the mitiaf conditions of cosmic chemicaf evofution, 

one should seek to understand the origin and relative abun- 
dances of the metals in the least chemically evolved systems. 

* Based on observations collected at the European Organisation 

for Astronomical Research in the Southern Hemisphere, Chile The most metal-poor damped Lya systems (DLAs), for 

[VLT programs 67A-0078(A), 69.A-0613(A), 083.A-0042(A), example, are usually interpreted as dist ant protogalaxies 

085.A-0109(A)], and at the W.M. Keck Observatory, which is op- at an early stage of chemical evolution (|Erni et all 120061 : 

erated as a scientific partnership among the California Institute of |Cooke et alj|201ll ). Whilst the origin of their metals is still 

Technology, the University of California and the National Aero- unknowrli recent hydrodynamical simulations sug- 

nautics and Space Administration. The Observatory was made , , , , , . , , , , . , , , 

., , , „ . , . T , , gest that such systems might have been enriched by just a 

possible by the generous financial support ol the W.M. Keck " j — -i i _ i " 

Foundation. Keck telescope time was granted by NOAO, through few supernova events ( |Bland-Hawthorn et al. | [201l|) . If this 

the Telescope System Instrumentation Program (TSIP). TSIP is is indeed the case > the most metal-poor DLAs provide a sim- 

funded by NSF. pl e route to study the first stages of chemical enrichment in 

f email: rcooke@ast.cam.ac.uk our Universe. 
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By definition, DLAs have a neutral hydr ogen column 
densi ty in excess of 2 x 10 20 Hi atoms cm l|Wolfe et al. 
19861 : see also the review bv IWolfe. Gawiser. fc Prochaska 



20051 ) . which acts to self-shield the gas from the ultravi- 



olet background radiatio n of quasars (QSOs) and galaxies 
(|Haardt fc Madaull200 lT ). This results in the gas having a 
simple ionization stru cture subject to neg ligible corrections 
for unseen ion stages (|Vladilo et al.ll200ll ). quite unlike the 
Lya fores t clouds that trace the low densi ty regions of the 
IGM (e.g. lSimcoe, Sargent, fc Raucrj|20blh . The main con- 
cerns that limit abundance studies in DLAs are line satu- 
ration and the possibility th at dust may hide some fraction 
of the metals (jVladild 12004 ). These concerns are alleviated 
when the metallicity of the DLA is below ~ 10 -2 Zq, which 
is also the regime where we expect to uncover the enrichment 
signature of the earliest generations of stars. 

The recent interest in the most metal-poor DLAs 
jPettini et alj|2008l : iPenprase et~alll2010l ) complements the 
ongoing local s tudies of metal-poor stars in the halo of 
the Milky Way dCavrel et alj|2004 l Beers fc Christ liebll200l 
ISudaet alj|2008l : lFrebej|201Cft . These stars are believed to 
have condensed out of near-pristine gas (perhaps a metal- 
poor DLA itself?), that was enriched by only a few earlier 
generations of stars. Thus, the first generation of stars can 
also be studied through the signature retained in the stellar 
atmospheres of the most metal-poor stars in the halo of our 
Galaxy. However, unlike the relative ease with which one 
can measure the abundances of metal-poor DLAs, deriving 
element abundances from the st ellar atmosphe res of metal- 
poor stars is not straightforward ((Asplund 200l)- Systematic 
uncertainties in the derived abundances are introduced by 
assuming that the spectral line being examined forms in a 
region that is in local thermodynamic equilibrium (LTE), 
as well as the need to account for three-dimensional (3D) 
effects in the ID stellar atmosphere models. 

These effects are particularly acute for oxygen, where 
several different abundance indicators are known to pro- 
duce contradictory estim ates in the low-metallicity regime 
|Garcfa Perez et al.ll2006T ). Despite the efforts of many au- 
thors, our uncertainty in the derived oxygen abundances has 
sparked an ongoing debate as to the trend of [O/FeQ in the 
Milky Way when [Fe/H] ,$ — 1.0. A history of the re levant 
discussion on [O/Fe] is provided bv lMcWilliaml (|l997t) , with 
further details given in Section [6.21 In brief, at low metal- 
licity, both O and Fe are produced exclusively by type-II 
supernovae (SNe II) and the winds from their progenitors. 
When [Fe/H] > - 1.0, there is a drop in [O/Fe] due to the 
delayed contribution of Fe from type-la supernovae (SNe la). 
Thus, the [O/Fe] ratio is most commonly used to measure 
the time delay between SNe II and the onset of SNe la. At 
the lowest metallicity, however, one can use the [O/Fe] ra- 
tio as a measure of the relative production of a- to Fe-peak 
elements by the first few generations of massive stars. 

Another key diagnostic ratio at low metallicity that may 
shed light on the nature of the early ge nerations of stars 
was uncovered by lAkerman et al.l |2004 ) who reported a 
rather surprising evolution of [C/O] with decreasing O abun- 



dance in their sample of 34 halo stars (see also ISpite et al.l 
2005). In disc and halo stars when the oxygen abundance is 
> - 1.0, [C/O] steadily rises from [C/O] ~ -0.5 to solar. 
When [O/H] ^ — 1.0, galactic chemical evolution models 
that only consider the nucleosynthetic products of Popu- 
lation II stars predict [C/O] to decrease or plateau, con- 
trary to the observed trend. The increase in [C/O] with 
decreasing metallicity has thus been interpreted as evi- 
dence for an increased carbon yie l d from either Population 
HI stars dChieffi fc Limongill2002l ; lUmeda fc Nomotd 120031 ; 
iHeger fc Wooslevll2010ri or rapidly-rotating low-metallicity 
Population II stars i Chiappini et al.l I2006L At first, con- 
cerns were raised regarding the accuracy of the derived C 
and O abundances, si nce the lines u s ed are subject to large 
non-LTE corrections. iFabbian etail l|2009aT ). however, per- 
formed a non-LTE analysis of the same lines, with further 
contraints from additional C I lines, to confirm the reality of 
the stellar [C/O] trend. These results depend somewhat on 
the adopted cross sections for collisions of C I and O I atoms 
with electrons and hydrogen atoms, but for all probable val- 
ues, [C/O] increases with decreasing metallicity when [O/H] 
< -2.0. 

To summarize, at present there are still some remaining 
concerns that prevent us from accurately measuring C and 
O abundances in the atmospheres of metal-poor halo stars. 
These difficulties have prompted a few teams to focus on 
very metal-pooi0 (VMP) DLAs where the absorption lines 
of Cn and Ol may be unsaturated and the abundances of 
C and O can be measured with confidence. Unfortunately, 
these near-pristine DLAs are rare, falling i n the tail of the 
meta llicity distribution function of DLAs (|Prochaska et al] 
I2OO7I) . Thus, only a handful of confirmed VMP DLAs are 
known at present. The first high spectral resolution survey 
(R ~ 40000, full width at half maxi mum, F WHM c± 7 km 
s _1 ) for VMP DLAs was conducted bv lPettini et alj i|200ct) . 
whose specific goal was to study the relative abundances of 
the CNO group of elements as a probe of early nucleosynthe- 
sis. Indeed, this was the first study to independently confirm 
the increased [C/O] abundance at low metallicity, suggest- 
ing that near-solar values of [C/O] are commonplace in this 
metallicity regime. 

The [C/O] trend reported by iPettini et al.1 d2008f) has 
also been independently noted bv IPenprase et al.1 (2010) m 
a medium spectral resolution (R ~ 5000, FWHM ~ 60 km 
s _1 ) survey of 35 DLAs (a prelimina ry report of this study 
can be found in IPenprase et aUl200&t ). In many of their sys- 
tems, the Cn and Ol lines were thought to be affected by 
line saturation, leaving only five DLAs to test the trend in 
C/O. Interestingly, this sample of DLAs suggests that [C/O] 
continues to rise to supersolar values when [O/H] ,$ — 3. 

Such surveys for VMP DLAs are most useful for study- 
ing the general properties of entire clouds of near-pristine 
gas before they form stars. In this contribution, we build on 
our ongoing survey for the most metal-poor D L As as probes 
of ea rly nucleosynthesis (|Pettini et alj 120081 ; ICooke et al.l 
l201ll ). With additional systems drawn from the literature, 
the total sample presented herein amounts to 22 VMP DLAs 



1 We adopt the standard notation: [A/B] = log^A/We) — 
log(A r A/A r B)© 1 where A^a.b refers to the number of atoms in ele- 
ment A and B. 



2 Herein, we adopt "very metal-poor" to be those DLAs with 
[Fe/H] < — 2.0, in line with the clas sification scheme for stars 
proposed bv lBeers fc Christliebl J2005h . 
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Table 1. Journal of Observations 



QSO 




2cm 


■^abs 


Telescope/ 


Wavelength 


Resolution 


Integration 


— 

S/N b 


Program ID 




(mag) 






Instrument 


Range (A) 


(km s — ) 


Time (s) 






J0311-1722 


17.7 


4.039 


3.73400 


VLT/UVES 


4370-6410 c 


6.9 


3 600 


15 


69.A-0613(A) d 


J0831+3358 


19.5 


2.427 


2.30364 


KECK/HIRESb 


3130-5970 c 


7.3 


18 600 


15 


A185Hb 


J1001+0343 


19.2 


3.198 


3.07841 


VLT/UVES 


3740-6650 c 


7.3 


33 700 


40 


083.A-0042(A) 


J1037+0139 


19.4 


3.059 


2.70487 


VLT/UVES 


3640-6650° 


7.3 


26 075 


45 


083A-0042(A) 


J1340+1106 


19.0 


2.914 


2.50792 


VLT/UVES & 


3483-6652° 


7.3 


29 800 & 


50 


085A-0109(A) & 








2.79583 


VLT/UVES & 


3483-9396 c 


10.3 


10 800 & 


10 


67.A-0078(A) d & 










KECK/HIRES 


4648-7044 c 


8.1 


8100 


25 


UllH d 


J1419+0829 


18.9 


3.030 


3.04973 


VLT/UVES 


3710-6652 c 


7.3 


29 800 


43 


085A-0109(A) 



a Magnitudes are SDSS g-band, except for J0311-1722 (not covered by the SDSS) which is ij-band llPeroux et alj|200ll) . 
b Indicative signal-to-noise ratio at 5000 A (or 6000 A in the case of J1340+1106). 
c With some wavelength gaps. 

d Spectra downloaded from either the UVES or HIRES data archives 4,5 . 



with abundance measurements derived from high spectral 
resolution data. From this sample, we confirm the elevated 
[C/O] values in these systems and, for the first time, present 
the trend of [O/Fe] in the most metal-poor DLAs. For both 
of these diagnostic ratios, we comment on the implications 
our findings have for local studies of Galactic metal-poor 
halo stars. Finally, we construct the abundance pattern of a 
typical VMP DLA for the elements C, N, O, Al, Si and Fe, 
and compare it to model calculations of Population II and 
Population III nucleosynthesis. 

This paper is arranged as follows. In Section[2]we detail 
the processing and preparation of the data. In Section [3] we 
explain the profile fitting procedure used for our new sample 
of VMP DLAs, which are discussed in Section [4] The accu- 
racy of our abundance analysis is discussed in Section [5l 
before we investigate the behaviour of [C/O] and [O/Fe] in 
VMP DLAs and compare with stellar data, in Section [6] 
Finally, we discuss the implications for our findings in Sec- 
tion [7] before summarizing our main results and drawing 
our conclusions in Section [8] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Target Selection 

Even at the relatively low spectral resolution afforded by 
the Sloan Digital Sky Survey (SDSS), one can easily recog- 
nise DLAs in the spectra of quasars, owing to the charac- 
teristic damping wings of the Lya absorption line profile. 
Subsequent identification of associated metal line absorption 
leads to a rough estimate of the gas-phase metallicity. Can- 
didate metal-poor DLAs are then identified as those DLAs 
that appear to exhibit no metal line absorption; their ab- 
sorption features are unresolved at the spectral resolution of 
the SDSS. However, when these candidates are re-observed 
with echelle spectrographs of high resolution (R Si 30000, 
FWHM ^ 10 km s _1 ), the metal absorption lines are re- 
solved, and in many cases it i s possible to measur e elemental 
abundances with confidence |Pettini et alj|2008h . 

The most recent trawls through SDSS spectra of ~ 8000 
quasars with z em Si 2.2 has yi e lded a sample of ~ 100 
DLAs (|Noterdaeme et alj 120091 ; iProchaska fe Wolfd 120091) 
of which ~ 400 are classified as 'metal-poor' ( Penprase et al.l 



2010)0 From compilations such as these, we selected a 
handful of metal-poor DLA candidates that exhibit no dis- 
cernible metal-line absorption at the spectral resolution of 
the SDSS, and re-observed these with echelle spectrographs, 
giving higher priority to candidates with: (i) bright quasars, 
so as to efficiently obtain spectra with signal-to-noise ra- 
tios S/N Si 20 in the continuum; (ii) DLAs where the dif- 
ference between z a b s and z cm is minimized so that the ab- 
sorption lines of interest (e.g. OlA1302 and CnA1334) are 
not blended with unrelated Lya forest lines; (iii) quasars 
whose emission redshift is below z cm ^ 3.3, so there is an 
improved chance that other lines of interest (e.g. Ol A1039 
and CnA1036) are not blended with Lya forest lines; (iv) 
DLAs at the low end of the column density distribution func- 
tion - that are still DLAs - to ensure that even the strongest 
metal absorption lines are unsaturated, allowing us to mea- 
sure the metal ion column densities with confidence; and (v) 
quasars with more than one metal-poor DLA candidate in 
their spectra. 

Our survey to date consists of 12 DLAs with [Fe/H] ^ 
—2.0. Initial result s for four of these were published by 
iPettini et all (|2008l ). while a fifth DLA, showing a pro- 
nounced C enhance ment relative to F e, was the subject 
of a recent study by ICooke et alj (|201ll ). The observations 
and analysis of the remaining seven DLAs, including one 
from the European Southern Observatory's (ESO) Ultravi- 
olet and Visual Echelle Spectrograph (UVES) data archive, 
are presented here. To our own data, we add a collection 
of published abundance measurements in ten VMP DLAs, 
selected as described in Section [4.81 to assemble an overall 
sample of abundance measurements in 22 VMP DLAs, all 
obtained from high resolution spectra (R Si 30 000). 



2.2 Echelle spectroscopic follow-up 

In order to achieve the high signal-to-noise ratio and spec- 
tral resolution required for accurate DLA abundance mea- 
surements, we observed our prime candidates with echelle 



3 In this context, a DLA is classed as 'metal-poor' if it has fewer 
than three significantly (4<r) detected metal absorption lines at 



the spectral resolution of the SDSS. 
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spectrographs on 8 — 10 m class telescopes. Most of our 
candidates were o bserved with the UVES spectrograph 
ijDekker et alJl200d ). which is mounted on UT2 at the Very 
Large Telescope facility. An additional system was observed 
with the W. M. Keck O bservatory's High Resolution Echelle 
Spectrometer fHIRES. IVogt et al.l[l994h on the Keck I tele- 
scope. Table [T] lists details of the observations of seven VMP 
DLAs reported here for the first time. For J1340+1106, we 
include details of some additional data, of comparable spec- 
tral resolution to ours, retrieved from the ESC0 and Keck 
Observatorjjf] data archives (program IDs 67. A-0078(A) and 
U11H respectively). We have also retrieved UVES spectra 
of the quasar J031 1 — 1722 fr om the ESQ data a rchive (pro- 
gram ID 69.A-0613(A); see IPeroux et alj 120051 ). since the 
metal lines for the VMP DLA along this sightline have not 
been previously analysed. 

Our own UVES observations [program IDs 083. A- 
0042(A) & 085.A-0109(A)] employed a 1.2" wide slit, re- 
sulting in a spectral resolution R ~ 41 000 (velocity FWHM 
~ 7.3 km s _1 ) sampled with ~ 3 pixels. We used dichroic 
#1 to split the quasar light into the blue and red spectro- 
scopic arms containing the HER_5 filter and SHP700 filter 
respectively. The resulting central wavelength for each arm 
was 3900 A (blue) and 5640 A (red). Both the blue- and red- 
sensitive CCDs used 2x2 on-chip binning. For our HIRES 
observations (program ID A185Hb) we used the C5 decker 
(a 7.0 x 1.148 arcsec slit) which, with su b-arcsec seeing, gav e 
a spectral resolution of R ~ 41000 fcf. ICooke et a.1.1 1201 IT ) . 
also sampled with ~ 3 pixels. We employed the ultravio- 
let cross-disperser with no filters, and used 2x2 on-chip 
binning. 

2.3 Data Reduction 

We used the standard UVES data reduction pipelined pro- 
vided by ESO to reduce the UVES data. The UVES re- 
duction pipeline performs the usual steps relevant to echelle 
data reduction. The preliminary steps include bias subtrac- 
tion, flat fielding, and background subtraction. The echelle 
orders are then traced using a flat field frame taken with 
a pinhole decker, and ID spectra extracted. The data are 
wavelength calibrated with reference to a ThAr lamp. 

The HIRES data were reduced with the makee data 
reduction pipeline developed by Tom Barlow, makee per- 
forms the same reduction steps as outlined above, but a trace 
frame is not always readily available. When available, the or- 
ders were traced using a flat-field frame taken with a pinhole 
decker. Otherwise, the science exposure of the quasar itself 
was used when a satisfactory trace could be made. Failing 
this, a trace frame was generated with a suite of purpose- 
built python programs, using the science frame as a guide 
to trace the echelle orders. 

Following these initial reduction steps, for each object 
we combined the science exposures using the software pack- 
age uves_poplefQ, maintained by Michael Murphy. This 

4 http:/ /archive. eso.org/eso/eso_archive_main. html 

5 https:/ /www2. keck. hawaii.edu/koa/public/koa.php 

6 We used version 4.3.0, available from: 
htt p : / / www .eso.org/sci/ software /pipelines/ 
' uves_popler can be downloaded from 

http: / /astronomy.swin. edu.au/ ~mmurphy/UVES_popler 



Table 2. Absorption Components of Low Ion Transitions 



Component 




b 


Fraction a 


Number 




(km s ) 


AT(Sill) 


J0311 


-1722: DLA at z abs = 


3.73400, x 2 /dof= 


0.80 


1 


3.733862 ± 0.000017 


5.6 ± 1.2 


0.18 


2 


3.733998 ± 0.000007 


2.7 ± 0.8 


0.23 


3 


3.734035 ± 0.000028 


14.4 i 1.5 


0.20 


1 


3.734439 ± 0.000002 


4.6 ± 0.2 


0.39 


J0831+3358: DLA at z abs = 


2 30364 v 2 /Hof— 


1.20 


1 


2.303565 ± 0.000004 


4.4 i 0.2 


0.49 


2 


2.303720 ± 0.000004 


6.1 ± 0.3 


0.51 


J1001 


f 0343: DLA at z abs = 


S 07841 v 2 /dof— 


1.42 


1 


3.078413 ± 0.000002 


7.0 ± 0.1 


1.00 


J1037 


f0139: DLA at z abs = 


2.70487, x 2 /dof= 


1.08 


1 


2.704870 ± 0.000002 


5.9 ± 0.2 


1.00 


J1340 


f 1106: DLA at z abs = 


2.50792, x 2 /dof= 


1.38 


1 


2.507649 ± 0.000003 


2.0 ± 0.4 


0.19 


2 


2.507921 ± 0.000001 


5.8 ± 0.1 


0.81 


J1340+1106: DLA at z abs = 


2.79583, x 2 /dof= 


1.76 


1 


2.7955454 ± 0.0000018 9.2 ± 0.1 


0.21 


2 


2.7958272 ± 0.0000007 6.55 ± 0.05 


0.79 


J1419 


f0829: DLA at z abs = 


3.04973, x 2 /dof= 


1.45 


1 


3.049649 ± 0.000002 


3.5 ± 0.1 


0.43 


2 


3.049835 ± 0.000002 


6.4 ± 0.1 


0.57 



a Fraction of the total column density of Sill. 



software merges individual echelle orders, and maps the data 
onto a vacuum heliocentric wavelength scale. Finally, we nor- 
malized the data by dividing out the quasar continuum and 
emission lines. Using the approximate redshift derived from 
each DLA's SDSS discovery spectrum, we then prepared the 
final data for analysis by extracting a ±150 km s _1 window 
around the pixel with highest optical depth near all available 
absorption lines of interest. Finally, a further fine adjustment 
to the continuum was applied to these extracted portions of 
the spectra when necessary. 



3 PROFILE FITTING 

For DLAs with a metallicity below 1/100 Zq, the metal 
line absorption is typically concentr ated in only a few 
clouds of low vel o city dispersion llLedoux et al.l 120061 ; 
iMurphv et al. 1 120071 ; iProchaska et al.l 120081 ). By assuming 
that a Maxwellian distribution accurately describes the ve- 
locities of the dominant atoms within the neutral cloud, we 
can model a DLA's absorption lines by a Voigt profile. To 
this end, we employed the Voigt profile fitting software vpfit 
to derive the cloud parameters for all DLAs in our sampled 
vpfit uses a chi-squared minimization algorithm to si- 
multaneously fit multiple Voigt profiles to a set of absorption 
lines characterized by three free parameters: (1) the cloud's 
absorption redshift (z aba ); (2) the Doppler parameter of the 
absorbing gas (b in km s _1 ); and (3) the column density of 



VPFIT is available from http: / /www. ast.cam.ac.uk/ ~rfc/vpfit.html 
© 2011 RAS, MNRAS 000.014261 
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Table 3. Adopted metal line laboratory wavelengths and oscillator strengths 



Ion 


Wavelength 

(A) 


/ 


Ion 


Wavelength 

(A) 


/ 


Ion 


Wavelength 

(A) 


/ 


Cn 


1036.3367 


0.118 


Al ii 


1670.7886 


1.740 


Fc ii 


1063.1764 


0.0547 


C ii 


1 334 5323 


0.1278 


Al in 


1854 71829 


0.559 


Fe II 


1081 8748 


0.0126 


Cn* 


1335.6627 


0.01277 


Al in 


1862.79113 


0.278 


Fen 


1096.8769 


0.0327 


C ii* 


1 335 7077 


0.115 


Si ii 


989 8731 


0.171 


Fe II 


1125.4477 


0.0156 


N I 


1134.4149 


0.0278 


Sin 


1020.6989 


0.0168 


Fen 


1143.2260 


0.0192 


N I 


1 1 34 9803 


0.0416 


Sin 


1 1 90 41 58 


0.292 


Fe II 


1144 9379 


0830 


N i 


1199.5496 


0.1320 


Sin 


1193.2897 


0.582 


Fen 


1260.533 


0.0240 


N I 


1200.2233 


0869 


gin 


1260.4221 


1.18 


Fc II 


1 608 4509 


0.0577 


N 1 


1900 7098 


0.0432 


Sin 


1 304 3702 


0863 


Fe II 


1611.20034 


001 38 


Nil 


1083.9937 


0.111 


Sin 


1526.7070 


0.133 


Fen 


2344.21296 


0.1142 


Oi 


925.446 


0.000354 


Sin 


1808.01288 


0.00208 


Fen 


2374.46033 


0.0313 


Oi 


936.6295 


0.00365 


Sn 


1250.578 


0.00543 


Fen 


2382.76418 


0.320 


Oi 


948.6855 


0.00631 


So 


1253.805 


0.0109 


Ni ii 


1317.217 


0.057 


Oi 


976.4481 


0.00331 


So 


1259.5180 


0.0166 


Ni ii 


1370.132 


0.056 


Oi 


988.5778 


0.000553 


An 


1048.2199 


0.263 


Ni ii 


1454.842 


0.0323 


Oi 


988.6549 


0.0083 


An 


1066.6598 


0.0675 


Ni ii 


1709.6042 


0.0324 


Oi 


988.7734 


0.0465 


Cm 


2056.25693 


0.1030 


Ni ii 


1741.5531 


0.0427 


Oi 


1039.2304 


0.00907 


Cm 


2062.23610 


0.0759 


Ni ii 


1751.9157 


0.0277 


Oi 


1302.1685 


0.048 















the ion that gives rise to the absorption line. When it was ev- 
ident that the DLA metal absorption arises from more than 
one cloud component, we introduced additional components 
to reduce the \ 2 (such that the \ 2 divided by the number of 
degrees of freedom [dof] was close to 1.0), whilst maintaining 
realistic errors on the derived parameters (i.e. ^ 10% uncer- 
tainty on b and a redshift uncertainty less than the sampling 
size of ~ 2.5 km s _1 ). Throughout the fitting procedure we 
assumed that the dominant ions in Hi regions (e.g. Cn, Nl, 
Oi, Sin, Fen) are kinematically associated with the same 
gas. We therefore fixed the redshift and Doppler parameter 
of each absorption component to be the same for each of 
these ions. The resulting cloud model parameters, including 
the reduced \ 2 y are provided in Table [2] where, as a guide, 
the last column lists the fraction of the total column den- 
sity of Sin in each component. The total column densities 
of available ions in each DLA are collected in Table Al. Ta- 
ble [3] lists laboratory wavelengths and oscillator stre ngths of 
releva nt atomic transitions from th e compilation b y Morton 



Table 4. ION COLUMN DENSITIES OF THE DLA IN J0311 — 1722 AT 

•Jabs = 3.73400 



relevant atomic transitions irom tnc compilation oy iviortoi 
(2003), with subsequent updates bv lJenkins fc Tripd l|200rf ^ 



4 INDIVIDUAL OBJECTS 

In this section, we briefly comment on the properties of each 
new DLA analysed in this paper. 



4.1 J0311 1722: DLA at z abs = 3.73400 



The VMP DLA along 
3 h 



the line-of-sight to J0311-1722 
(J 2000.0: 03''ll' n 15!2q , -17° 22' 47'.' 4) was first identified 
by Peroux et all (|200lh . Follow-up LIVES spectroscopy by 
IPeroux et al. (2005) revealed a Lya absorber at z a b s = 3.734 
with log [iV(Hl)/cm" 2 ] = 19.48 ± 0.10, which is lower tha n 
the conventional limit for DLAs set by I Wolfe et all l|l98rj) . 
log [N(R i)/cm~ 2 1 ^ 20.3. Such syst ems are often referred 
to as sub-DLAs jPeroux et al. 2003a) or super Lyman-limit 



Ion 


Transitions used 


log A r (X)/cm- 2 


Hi 


1025, 1215 


20.30 ± 0.06 


Cii 


1036, 1334 


14.02 ± 0.08 


Cn* 


1335 


13.55 ± 0.06 


Ni 


1200.2 


< 13.07 a 


Oi 


1302 


14.70 ± 0.08 


Sin 


1260, 1304 


13.31 ± 0.07 


Fen 


1125 


< 13.76 a 



L 3<t limiting rest frame equivalent width. 



systems, and are defined to have 19.0 ^ log [iV(H l)/cm 2 ] ^ 
20.3 

IPeroux et al.l (120051 ) derived their estimate of AT(Hi) 
for this absorber from a consistent model fit to the first 6 
Lyman lines (from Lya to Ly6), where the Lya line produces 
the poorest fit. Having re-analysed these data, giving higher 
priority to the Lya and Ly/3 lines, we find that a better fit 
results if log [AT(Hl)/cm~ 2 ] = 20.30 ± 0.06 (see top panel of 
Figure [TJ, provided that the Doppler paramete r is 6 ^ 30 km 
s" 1 . P ossibly, the value of A'' (Hi) derived by IPeroux et al.l 
(2005) was driven by the higher order Lyman lines, which 
may be blended with other Lyman forest lines, and resulted 
in a Doppler parameter much larger than 30 km s" 1 . 

In the lower panels of Figure Q] we present a selection of 
the metal lines associated with this VMP DLA. All of the 
available metal absorption lines are unsaturated, leading to 
reliable estimates of the metal column densities. A four com- 
ponent model was found to accurately reproduce the metal- 
line profiles, whilst maintaining reasonable estimates of the 
parameter errors (see Section [3j . The derived cloud model 
parameters are listed in Table [5] Although both C II lines 
(A1334 and A1036) are blended on the red wing of compo- 
nent 4 (centred near v = km s _1 in Figure [TJ, this does not 
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Figure 1. The top panel shows the DLA towards J0311— 1722 (black histogram) which exhibits a damped Lyo line at z^ s = 3.73400. 
The red continuous line shows the theoretical Voigt profile for a neutral hydrogen column density log[7V(H l)/cm~ 2 ] = 20.30. The 
remaining panels display a selection of metal lines, overlaid with the Voigt profiles for the derived cloud model (in red) and blends (in 
light blue). The red tick marks above the normalized continuum indicate the locations of the absorption components. For all panels, the 
y-axis scale is residual intensity. The normalized quasar continuum and zero-level are shown by the blue dashed and green dotted lines 
respectively. 



greatly affect our final estimate of iV(Cll), since component 
4 has a well-determined Doppler parameter from the relative 
strengths of SillA1260 and SillA1304, and from other un- 
blended transitions in regions of high signal-to-noise. In any 
case, the majority of the absorbing column is contributed 
by the first three components (~ 60%). We list the total 
column density returned by vpfit for each available ion in 
Table H 

In this table, we also provide upper limits for the col- 
umn densities of several key ions that are undetected at the 
signal-to-noise of the data. Specifically, we calculate the 3<r 



limiting rest-frame equivalent width, Wo, over the velocity 
interval of absorption exhibited by the weakest transition, 
which in this case is SillA1304. A 3a upper limit to the un- 
detected feature is then derived using the optically thin limit 
approximation, N = 1.13 x 10 20 ■ Wo/A 2 / cm" 2 . For Nl we 
use the undetected A1200.2 line to derive W (Nl) < 13 mA, 
which implies log JV(Nl)/cm" 2 < 13.07. Similarly for 
FenA1125, Wo(Fell) < 10mA implies log A r (Fe ll)/cm -2 
13.76. 
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Figure 2. Same as Figure[T]for a selection of metal lines associated with the 2 a b s = 2.30364 DLA towards J0831+3358. 



Table 5. Ion column densities of the DLA in .10831+3358 at 
2 aba = 2.30364 



Ion 



Transitions used 



log 7V(X)/cm" 



Hi 1215 20.25± 0.15 a 

Ni 1199.5 < 12.78 b 

Oi 988.5,988.6,988.7,1039,1302 14.93 ± 0.05 

Alii 1670 12.19 ± 0.06 

Sin 989, 1193, 1260, 1304 13.75 ± 0.04 

Sn 1250 s: 13.75 b 

Fen 1608 13.33 ± 0.06 



a Penprase et al. (2010) 

b 3<r limiting rest frame equivalent width. 



4.2 J0831+3358: DLA at z abs = 2.30364 

We observed J0831+3358 with HIRES on 2009 December 9 
under good conditions with subarcseco nd seeing. We used a 
1.148" wide slit which, as measured bv lCooke etall <|201ll ). 



delivered a spectral resolution of 7.3 km s _1 FWHM. Un- 
fortunately, the Lya line at A b s — 4015A falls in a gap 
between two of the CCDs on the HIRES detector mo- 
saic. Thus, for this DLA, we adopt the Hi column density 
logJjV(Hl)/cm~ 2 ] = 20.25 ± 0.15 derived bv lPenprase et all 
(2010) from their observations of this QSO at R ~ 5000, 
which is sufficient to resolve the broad damped profile of 
the Lya line. 



The metal-lines in our data are well fit by a model with 
two components of roughly equal strength separated by 14 
km s~ x (see Figure[2]). Details of the derived cloud model are 
presented in Table [5] with the associated column densities 
given in Table [5] The C II lines are saturated in this DLA, 
however, we have a clean measurement of A^(0 1) from a 
number of unsaturated Oi lines. In Table [S] we also provide 
3a upper limits to the Ni and Sn column densities which 
are undetected at the signal-to-noise of our data. 
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Figure 3. Same as Figure[T] for the DLA towards J1001+0343 (black histogram) which exhibits a damped Lya line at z a b s = 3.07841 
(top panel). Here, the red continuous line shows the theoretical Voigt profile for an Hi column density log[iV(Hl)/cm" 2 ] = 20.21. The 
remaining panels display a selection of metal lines. Note the different y-axis scale that is used for the weak Fell A1608 line (bottom-right 
panel). 



Table 6. ION COLUMN DENSITIES OF THE DLA IN J1001+0343 AT 
Zabs = 3.07841 



Ion Transitions used log 7V(X)/cm 2 

Hi 1215 20.21 ± 0.05 

Cn 1036, 1334 13.58 ± 0.02 

Ni 1200.2 < 12.50 a 

Oi 1039, 1302 14.25 ± 0.02 

Sin 1190, 1193, 1260, 1304, 1526 12.86 ± 0.01 

Sn 1253 ^ 12.91 a 

Fen 1608 12.50 ± 0.14 



a 3cr limiting rest frame equivalent width. 



4.3 J1001+0343: DLA at z abs = 3.07841 

This QSO was observed with UVES in service mode on the 
nights of 2009 April 19 & 29, 2010 January 11 & 27 and 
2010 February 7. Our total integration time was 33700 s, 
yielding a S/N per pixel of ~ 40 at 5000 A The DLA in- 
line to J1001+0343 was also investigated bv lPenprase et all 
(2010), being amongst the most metal-poor in their sam- 
ple. From our observations we derive an Hi column density 
of log[A^(Hl)/cm~ 2 ] = 20.21 ± 0.05 from the wings of the 
Ly a absorption line. Th is compares well with the estimate 
bv lPenprase et all |2010l ). log[JV(Hi)/cm -2 ] = 20.15 ±0.10. 
We present our Voigt profile fit to the Lya line in the top 
panel of Figure [3] 
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Figure 4. Same as Figure[l] for the DLA towards J1037+0139 (black histogram) which exhibits a damped Lya line at z a b s = 2.70487 
(top panel). The red continuous line shows the theoretical Voigt profile for an Hi column density log[jV(H i)/cm~ 2 ] = 20.50. The 
remaining panels display a selection of metal lines. 



Table 7. Ion COLUMN DENSITIES OF THE DLA IN J1037+0139 AT 
2 abs = 2.70487 



Ion 


Transitions used 


log iV(X)/cm- 2 


Hi 


1025, 1215 


20.50 ± 0.08 


Ni 


1199.5, 1200.2, 1200.7 


13.27 ± 0.04 


Oi 


1039, 1302 


15.06 ± 0.04 


Aln 


1670 


12.32 ± 0.03 


Sin 


1260, 1304 


13.97 ± 0.03 


Fen 


1143, 1144, 1608 


13.53 ± 0.02 



The remaining panels of Figure [3] showcase a number of 
the available absorption lines that were used in deriving the 



cloud model. In fact, all of the available metal absorption 
lines are unsaturated in this DLA, thus providing reliable 
measurements of the elemental abundances. Several metal- 
line transitions of varying strength are well fit by a single 
component cloud model with a Doppler parameter of b — 
7.0 ± 0.1 km s -1 (s e e Tab le [H). The curve of growth analysis 
by IPenprase et al.l |2010i) yielded a Doppler parameter of 
7.5 km s -1 , which is in good agreement with that found 
here. The derived column densities for all available ions are 
presented in Table [6] We also provide 3<r upper limits to the 
N I and S II column densities which are undetected at the 
S/N of our data. FeiiA1608 is detected at the 3.6a level. 
The corresponding fit is presented in the bottom-right panel 
of Figure [3] (note the different t/-axis scale). 
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Figure 5. Same as Figurefl] for the DLA towards J1340+1106 (black histogram) which exhibits a damped Lya line at z a i, s = 2.50792 
(top panel). The red continuous line shows the theoretical Voigt profile for an Hi column density log[Af(Hi)/cm" 2 ] = 20.09. The 
remaining panels display a selection of metal lines. 



4.4 J1037+0139: DLA at z abs = 2.70487 



This QSO was also observed in service mode with UVES 
on 2010 February 12-14 and again on 2010 March 5. 
J1037+01 39 was also one of the QSOs observed indepen- 
dently by IPenprase et al.l |2010t ). It is one of the faintest 
QSOs in our sample, requiring a total of 26075 s of in- 
tegration to achieve a S/N~ 40 at 5000 A. The Lya 
line falls on the edge of the blue detector, but this does 
not affect the accuracy of the Hi column density since 
the blue wing of the damped Lya line is still intact, 
and can be fit using the redshift derived from the well- 
defined narrow metal absorption lines. Moreover, we have 
access to Ly/3, which also exhibits damping wings (al- 
though not as strong as that of Lya). We derive an 



Hi column density of log[iV(H l)/cm~ 2 ] = 20.50 ± 0.08, 
which is consistent with estimates deri ved from the SPSS 
spect rum (log[AT(H i)/cm" 2 ] = 20. 45; iProchaska fe Wolfe] 
2009), as well as that derived by IPenprase et al.l 1 2010h 
(log[JV(Hi)/cm" 2 ] = 20.40 ± 0.25). We present our Voigt 
profile fit to the Lya line in the top panel of Figure [4] 

Again, the metal absorption is concentrated in a single 
component, in this case with a Doppler parameter of 5.9 
km s _1 . The corresponding column densities for all of the 
available ions are listed in Table [7] For this system, both 
C II A1036 and A1334 are saturated and blended; however, we 
have a robust measure of the O I and Fe II column densities 
from several unsaturated transitions. We also have a clear 
detection of the N I triplet near Ao = 1200 A. Selected metal 
absorption lines are reproduced in Figure [4] 
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Table 8. ION COLUMN DENSITIES OF THE DLA IN J1340+1106 AT 

z abs = 2.50792 



Ion 


Transitions used 


log 7V(X)/cm- 2 


Hi 


1215 


20.09 ± 0.05 


Ni 


1199.5, 1200.2 


12.80 ± 0.04 


Oi 


1039, 1302 


15.02 ± 0.03 


Alii 


1670 


12.27 ± 0.02 


Aim 


1854, 1862 


11.51 ± 0.10 


Sin 


1190, 1193, 1304, 1526 


13.75 ± 0.02 


Fen 


1096, 1143, 2344, 2374, 2382 


13.49 ± 0.02 



4.5 J1340+1106: DLA at z abs = 2.50792 

This QSO has previously been observed with both 
UVES (at a spectral resolution of 10.3 km s" 1 FWHM; 
iLedoux. Petitiean. fc Srianandir2003l) and HIRES (at a spec - 
tral resolution of 8.1 km s" 1 FWHM; [Prochaska et ai]|2003h . 
However, given that this QSO intersects two VMP DLAs, 
one of which had potentially unsaturated Cn lines, we de- 
cided to reobserve it with UVES for 29800 s at a slightly 
higher spectral resolution of 7.3 km s , and obtained com- 
plete spectral coverage from 3500 A to almost 1/Ltm. 

Since the broad damped profile of the Lya line is inde- 
pendent of the spectral resolution, we combined the three 
datasets to obtain a high S/N ratio near the damped Lya 
line, from which we derived log [AT (H i) /cm" 2 ] = 20.09±0.05. 
This model fit, along with the combined data, is shown in 
the top panel of Figure [5] The profiles of the metal absorp- 
tion lines, on the other hand, are narrow (FWHM ,$ 10 km 
s" 1 ); therefore, their observed profiles are not independent 
of the spectral resolution. Thus, we separately combined the 
data of equal spectral resolution and individually read these 
three reduced spectra into vpfit, which convolved the fit- 
ted model with the spectral resolution appropriate to the 
data. The upshot of proceeding in this way is that the cloud 
model is then largely driven by the dataset of highest S/N 
for each absorption line that is input. However, since we can- 
not combine all three datasets, in Figure [5] we only present 
the dataset with the highest S/N near each absorption line. 

Unfortunately, both of the C II absorption lines are sat- 
urated, and partially arise from nearby ionized gas, exhibit- 
ing a similar profile shape to the AlllA1670 line. In fact, 
most of the metal absorption lines exhibit a second com- 
ponent on the blue wing (at v = —23 km s" 1 relative to 
the main component at z abs — 2.507921), which appears to 
arise from ionized gas. This is confirmed by the absence of 
this blue component in O I absorption (see O I A1302, and 
note that the absorption on the blue wing of the OlA1039 
line is due to unrelated absorption), wh ich has long been 
know n to accurately trace neutral gas dField fc Steigmar] 
1971). The presence of ionized gas is also confirmed by the 
higher Allll/Alu ratio exhibited by the blue component. 

We derived the cloud model for both of these compo- 
nents from the host of available Si II and Fe II lines, with ad- 
ditional constraints coming from the two O I lines. In Tabled 
we present the fitting results from both the blue component 
(component 1; which arises from nearby mildly ionized gas), 
and the main component (component 2) which we attribute 
to the DLA. Fixing the parameters of this cloud model, we 
then derived the column density for all ions in both compo- 



Table 9. ION COLUMN DENSITIES OF THE DLA IN J1340 + H06 AT 



^abs 


2.79583 




Ion 


Transitions used 


log Af(X)/cm" 2 


Hi 


1215 


21.00 ± 0.06 


Cn* 


1335 


12.99 ± 0.04 


N I 


1134.4, 1134.9, 1199.5, 1200.7 


14.04 ± 0.02 


N ii 


1083 


12.81 ± 0.09 


Oi 


925, 976, 1302 


16.04 ± 0.04 


Alii 


1670 


13.24 ± 0.03 


Aim 


1854, 1862 


12.19 ± 0.08 


Sin 


1020, 1193, 1260, 1304, 1526, 1808 


14.68 ± 0.02 


Sn 


1253, 1259 


14.30 ± 0.02 


An 


1048, 1066 


13.18 ± 0.02 


Cm 


2056, 2062 


12.62 ± 0.11 


Fen 


1063.1, 1081, 1096, 1125, 1608, 


14.32 ± 0.01 




1611, 2260, 2344, 2382 




Nin 


1317, 1370, 1454, 1709, 1741, 1751 


13.08 ± 0.03 



nents. In Table[8] however, we only provide the column den- 
sity for the single, dominant component that we attribute 
to the DLA. The model fits to the data are presented in 
the lower panels of Figure [5] where, as stated above, we only 
present the model and data that correspond to the highest 
S/N for each absorption line. 

The detection of ions that arise in ionized gas, such as 
Aim, will later provide a useful means to test the accuracy 
of one of our underlying assumptions; that we can use the 
single dominant ion for each element to measure the elemen- 
tal abundances in DLAs (see Section 15.10 . 

4.6 J1340+1106: DLA at z abs = 2.79583 

We now report on the second VMP DLA that is inter- 
sected by this QSO. Again, we treat the Lya line and the 
metal lines of this DLA as detailed in Section 14.51 From 
the combined data we derive an Hi column density of 
log[iV(Hi)/cm~ 2 ] = 21.00L0.06. The combined spectrum in 
the region of the damped Lya line, together with the profile 
fit, is reproduced in the top panel of Figure [6] 

Turning now to the metal absorption lines, we have 
found that a two component cloud model (with Doppler pa- 
rameters of 9.2±0.1 and 6.55L0.05 km s _1 separated by 22 
km s _1 ) provides a good fit to the data. This cloud model 
is perhaps the best determined in our dataset, given the 
high S/N of the data, and the numerous atomic transitions 
available. Our Voigt profile fits to the metal lines are shown 
in the lower panels of Figure |S] However, as discussed in 
Section \4. 51 we only present the model and data that corre- 
spond to the highest S/N for each absorption line. The col- 
umn densities for all available ions are provided in Table [5] 
In this table, we have also provided the column densities for 
Aim and Nil, which are coincident with the DLA, but are 
typically associated with Hn regions. 



4.7 J1419+0829: DLA at z abs = 3.04973 

We recorded the spectrum of J1419+0829 (z cm = 3.034) 
for 29800 s with UVES in service mode, resulting in a S/N 
near 5000A (near the red wing of the damped Lya line) of 
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Figure 6. Same as Figure[T] for the DLA towards J1340+1106 (black histogram) which exhibits a damped Lya line at z a i, s = 2.79583 
(top panel). The red continuous line shows the theoretical Voigt profile for an Hi column density log[Af(Hi)/cm" 2 ] = 21.00. The 
remaining panels display a selection of metal lines. 



~ 40. This high S/N in combination with a virtually uninter- 
rupted red wing to the Lya absorption allows a very accurate 
measurement of the H I column density, log[iV(H l)/cm -2 ] = 
20.40 ±0.03. The Voigt profile fit to the Lya line is shown in 
the top panel of Figure [7] with a selection of the associated 
metal absorption lines in the remaining panels. 



A cloud model with two components separated by ~ 14 
km s _1 (Table [2} provides a good fit to the data. This cloud 
model is well determined by a host of O I and Si II lines (see 
Table [ID). Both Cn A1334 and A1036 are saturated; column 
densities for Ni, Ol Sin, and Fell are listed in Table [TOl 



4.8 The final VMP DLA sample 

In order to augment our survey of VMP DLAs, we have 
searched the literature for known examples satisfying the fol- 
lowing conditions: (i) the QSO spectra were observed at high 
spectral resolution (R > 30 000) - in practice this meant 
that the data were recorded with either UVES or HIRES; 
(ii) [Fe/H] ^ —2.0; and (iii) at least one unsaturated Ol ab- 
sorption from which [O/H] could be measured. These con- 
ditions were imposed to select a sample of measurements 
from the literature which is highly compatible to our own 
data and whose metal abundances could be adopted without 
reanalysing the spectra (although we referred all measure- 
ments to the same solar abundance scale - see Section [5). 
The literature trawl yielded an additional ten DLAs satisfy- 
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Figure 7. Same as Figure[l] for the DLA towards J1419+0829 (black histogram) which exhibits a damped Lya line at z^ s = 3.04973 
(top panel). The red continuous line shows the theoretical Voigt profile for an Hi column density log[Af(H l)/cm~ 2 ] = 20.40. The 
remaining panels display a selection of metal lines. 



ing the above conditions; together with our own observations 
they form a sample of 22 VMP DLAs. The metallicity dis- 
tribution function for this sample, which appears to tail off 
towards the lowest metallicities, is shown in Figure [8] 

Relevant details of the full sample are collected in Ta- 
ble [Tl] where we list absorption redshifts, neutral hydrogen 
column densities and element abundances for a selection of 
the metals most commonly observed in VMP DLAs, includ- 
ing: (1) C, N, O - the first elements synthesised in the chain 
of stellar nucleosynthesis; (2) Al - an odd atomic number 
element; (3) Si - an even atomic number element; and (4) 
Fe - an iron-peak element. The uncertainty in each abun- 
dance includes the error in Hi. In Appendix^] we provide 



Table 10. ION COLUMN DENSITIES OF THE DLA IN J1419+0829 

at z abs = 3.04973 



ion 


Transitions used 


log Af(X)/cm" 2 


Hi 


1025, 1215 


20.40 ± 0.03 


Ni 


1199.5, 1200.2, 1200.7 


13.28 ± 0.02 


Oi 


936, 948, 976, 988.5, 988.6, 988.7, 


15.17 ± 0.02 




1039, 1302 




Sin 


989, 1190, 1193, 1260, 1304, 1526 


13.83 ± 0.01 


Fen 


1063.1, 1608 


13.54 ± 0.03 



a similar table listing the column densities of each ion from 
which these abundances were derived. 
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Table 11. C, N, O, Al, Si, and Fe Abundance Measurements in VMP DLAs 
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Ref. a 


Our- A/A/TT} 
W 111 V IVlJr 


DLA Sample 
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lEllison et al. N2OIOI') quote a 3cr upper limit to the Fell column density of log JV(Fe 11) /cm 2 < 12.73. We have since rereduced these 
data (as described in Section [2J, and detected the FenA1608 line at the Act level. The Fe abundance quoted here is derived using the 
optically thin limit approximation to measure TV(Fell). 

c The metal ion uncertainty for this measurement was not provided by the authors. Thus, we only quote the uncertainty in A^Hl). 
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Figure 8. The metallicity distribution function for the sample of 
VMP DLAs listed in Table [TT] is shown as the red histogram. 

5 ABUNDANCE ANALYSIS 

As in previous DLA work, we assume that each element re- 
sides in a single dominant ionization stage in the neutral 
gas. Thus, the abundance of a given element is found by 
taking the ratio of the dominant ions column density to 



Table 12. Adopted solar abundances' 1 



X 


log(X/H) 


X 


log(X/H) 


c 


-3.57 


s 


-4.86 


N 


-4.17 


Ar 


-5.60 


O 


-3.31 


Cr 


-6.36 


Al 


-5.56 


Fe 


-4.53 


Si 


-4.49 


Ni 


-5.79 



a From Asplund et al. (2009). 



that of Hi, and referring it to a solar scale (i.e. [X/H] = 
log(Af(X)/Af( Hl)) - log (X/H)^). Throughout this article, 
we adopt the lAsplund et ail l|2009l 'l solar scale, taking the 
photospheric, meteori tic, or the average of th e two, based 
on the suggestion by lLodders. Plame. fc Gaill {2009) . The 
adopted abundances relevant to this work are collected in 
Table IT21 

Thus, measuring elemental abundances in VMP DLAs 
is a relatively straightforward process. Some uncertainty 
arises, however, if the observed metal-line absorption from 
ions that are dominant in H I regions, does not perfectly 
trace the Hi gas. For example, if the dominant ion for a 
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Figure 9. Ionization corrections for two typical VMP DLAs, respectively with low and high Hi column density. Left Panels: The 
ionization corrections, as defined by Equation [l] for the most commonly observed elements in VMP DLAs are plotted as a function of 
gas volume density. The vertical dashed lines correspond to the estimated gas densities of the two DLAs (shown by the solid vertical 
line in the right panels). Right Panels: Column density ratios for successive ion stages of Al (blue squares connected by a solid line) and 
N (red diamonds connected by a solid line) vs. gas density. We also plot the observed values of the Al ni/Al II and Nn/Nl ratios (solid 
horizontal lines) along with their uncertainties (dotted lines) for the two VMP DLAs we use as examples in this test (see text for further 
details). 



given element is also present in nearby Hn gas (often the 
case for singly ionized species; e.g. Cn, Sin, Fell), we will 
over-estimate the abundance of this element. On the other 
hand, if the dominant ion for a given element is mildly ion- 
ized in the H I gas itself (usually the case for neutral species; 
e.g. N i), we will under-estimate the element's abundance. In 
addition to these ionization corrections, further uncertain- 
ties may be introduced into the abundance analysis if some 
refractory elements have condensed to form dust grains. 
Whilst both of these co ncerns are expected to be neglig i- 
ble in the VMP regime jVladilo et al.ll200ll ; IVladild [2004) , 
we reassess their importance in the following sub-sections. 

5.1 Ionization Corrections 

Ionization corrections are known to be small when the neu- 
tral hydrogen column density is in excess of ~ 10 20 atoms 
cm~ 2 . Nevertheless, for each element X with the dominant 
ionization stage N, one needs to apply a small correction, 
IC(X), to recover the true elemental abundance, 

[X/H] = [Xn/Hi] +IC(X). (1) 

To estimate the magnitude of such corrections, we 
used the CLOUDY photoionization software developed by 



iFerland et all l| 19981 ) to model two of the VMP DLAs from 
our sample, chosen to represent the range of A^(H i) values 
that we report: one DLA with a low Hi column density (the 
DLA with log[Af(Hl)/cm" 2 ] = 20.09 towards J1340+1106 
at z abs = 2.50792) and the other with a high A r (Hl) (the 
DLA with log[Ar(Hl)/cm- 2 ] = 21.00 towards J1340+1106 
at z a bs = 2.79583). In both cases, we modelled the VMP 
DLA as a plane-parallel slab of constant volume density 
gas in the range —3 < log[n(H)/cm -3 ] < 3, irradiated 
by the cosmic microwav e background and UV background 
IIHaardt fc Madau|[200ll ) at the appropriate redshift. Using 
the solar abundance scale in Table [T2l we globally scaled the 
metal abundances of the VMP DLA to be 10~ 2 Z Q . The 
simulations were stopped once the H I column density of the 
DLA was reached, at which point we output the simulated 
ion column densities of the slab. 

Once the above value of the background radiation field 
is assumed, the ionization correction for each element de- 
pends on the volume density of the gas (left panels of Fig- 
ure [9}, which may be estimated by considering the ratio of 
successive ion stages (right panels of Figure [9]). For the low 
A'(Hl) DLA being considered as an example here, we mea- 
sure N(Aliii)/iV(Alii) = -0.76 ± 0.10, implying a gas den- 
sity of log[n(H)/cm" 3 ] ~ -1.2 ± 0.4. For the high N(Ul) 
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DLA, we measure JV(Alili)/jV(Aln)= -1.05 ± 0.09 and 
AT(Nii)/AT(Ni)= -1.23 ± 0.09, which are both consistent 
with a gas density of log[n(H)/cm~ 3 ] ~ —1.1 ±0.2. At these 
values of the gas density, it can be seen that the ionization 
corrections for the VMP DLAs in our sample are 0.1 dex 
for the main elements of interest. 

We also performed the above calculations under the 
assumption of a softer ba ckground radiation field (i.e. an 
O star) rather than the lHaardt fe Madaul (|200lf ) back- 
ground. Our results are quantitatively similar to those of 
IVladilo et all (|200lf ): corrections for all of the elements we 
are interested in are ^ 0.1 dex, except for Alii which can 
require corrections of the order ~ 0.3 dex for the lowest col- 
umn density systems. At present, it is not yet clear whether 
these very metal-poor systems harbour (or are nearby to) 
massive Population II or Population III stars. Based on the 
results described in this subsection, we have therefore not 
corrected any of the measured abundances for ionization ef- 
fects. 



5.2 Dust Depletion 

To measure accurately element abundances in DLAs, the 
fraction of a given element that is not observed in the gas 
phase, but is instead locked up in dust grains, must also 
be considered. To account for this effect, one ideally con- 
siders the relative abundances of a refractory and a volatile 
element (e.g. [Cr/Zn]), and compares this to the expected 
intrinsic nucleosynthetic ratio (typically the ratio seen in 
stars of comparable metallicity) . Previous studies based on 
such a comparison have shown that DL As exhibit m i nimal 
dust depletion when [Fe/H] ^ - 2.0 (|Pettini et all 1 19971 ; 
lAkermane t al. 2005). Unfortunately, the Cm and Znll lines 
are too weak in the VMP regime to be measured, so one 
needs to resort to more abundant elements, such as Si and 
Fe, which are known to be depleted to different degrees (Fe 
is more readily incorporated into dust grains than Si). 

The most metal-poor stars in the halo of our Galaxy 
suggest that the intrinsic nucleosynthetic ratio of Si/Fe is 
virtually independent of metallicity, co rresponding roughl y 
to a constant of [Si/Fe] = +0.37 ±0.15 (ICavrel et al.|[2 004). 
Such a plateau was first seen in DLAs bv lProchaska fe W olfe 
l|2002i ). who found [Si/Fe] ~ ±0.3 when [Fe/H] < -2 (see 
updated version in IWolfe. Gawiser fc Prochaskal 12001 ) . A 
similar study was also conducted bv l Vladilol (|2002l ') . who sug- 
gested there may still exist some mild depletion onto dust, 
resulting in a plateau of [Si/Fe] ~ ±0.25. From the 19 VMP 
DLAs in our sample, we find that [(Si/Fe)] = ±0.32 ± 0.09, 
which is certainly consistent with minimal dust depletion. 
We therefore proceed under the assumption that dust has a 
negligible effect on our derived elemental abundances. 



6 COMPARING VMP DLAS AND STARS 

DLAs likely experience very different chemical histories to 
that of the stars in the halo of our Galaxy. However, in 
the limit of decreasing metallicity, both stars and DLAs 
are polluted by very few previous generations of stars. In 
fact, since the physical conditions of the gas that gives rise 
to DLAs are conducive to fo rming stars (|Noterdaeme et all 
120081 : IJorgenson et al.l 120091 ). it certainly seems plausible 



that some of the VMP stars in the halo of our Galaxy origi- 
nally condensed out of a VMP DLA. If this is indeed true, we 
would therefore expect both VMP stars and DLAs to share 
similar chemical signatures at the lowest metallicities. In 
the following subsections, we test the validity of this expec- 
tation, by directly comparing the abundances of the VMP 
DLAs in this survey to the most recent abundance measure- 
ments of VMP stars. 



6.1 Revisiting C/O at low metallicity 

We first consider the trend of C/O at low metallicity, which 
has received a great deal of attention in recent years. The 
evolution of the C/O ratio when [O/H] £ — 1.0 has been 
known for a while; [C/O] increases linearly from ~ —0.6 to 
solar with increasing [O/H] (see Figure [TO)) . This trend is 
thought to be due to the increased, metallicity-dependent, 
carbon yields of massive rotating stars, combined with the 
delaye d release of carbon fr om low and intermediate mass 
stars jAkerman et al.l 120041 '). Based on current models of 
Population II nucleosynthesis, below [O/H] ~ —1.0, [C/O] is 
predicted to decrease (or perhaps plateau) with decreasing 
metallicity. Indeed, su ch a plateau was firs t reported for a 
sample of halo stars bv lTomkin et al.l |l992), who measured 
[C/O] abundances down to [O/H] ~ -1.7. 

This work was extended to even lower oxyge n abun- 
dances by lAkerman" et all (|2004T ) and ISpite et all (|2005h . 
Contrary to the supposed decrease in [C/O], these authors 
uncovered quite the opposite trend when [O /H] £ — 2.0; an 
extrapolation of this trend suggests that [C/O] could reach 
near-solar values when [O/H] ~ —3.0. Three possibilities 
have been suggested to explain this behaviour: (1) the left- 
over signature of a high-carbon producing generation of Pop- 
ulation HI stars (IChieffi fc Limongill2002l ; lUmeda fc Nomotol 



l2003l;lHeger fc Wooslevll2010l ); (2) pollution from a previous 
gener ation of rapidly-rotati ng low-metallicity Population II 
stars (|Chiappini et alj|2006r i; or (3) systematic uncertainties 
in the adopted ID LTE analysis. Th e third of these possibil- 
ities has recently been ruled out bv lFabbian et all (|2009al ). 
wh o conducted a more d etailed analysis of the lines used 
by lAkerman et all ([2004), thus confirming the stellar C/O 
trend at low metallicity. 

Such an 'unexpected' trend is perhaps not so surpris- 
ing in hindsight, since several other lines of evidence sup- 
port a high-carbon producing generation of early stars, 
including: (1) the fact that a high carbon abundance is 
required at early times to efficiently cool the gas, and 
drive the transi tion from Population HI to Pop ulation II 
star formation (|Frebel Johnson, fc Brommll2007h : (2) ob- 
servations of the three most iron-poor stars known to 
date have reveale d that they all exhibit extreme car- 
bon enhancement s (|Christlieb et al.ll2002l ; lFrebel et alll2005l ; 
iNorris et al.ll2007l) ; (3) in addition, the fraction of all carbon- 
enhanced metal-poor (CEM P) stars is thought to inc rease 
with decreasing metallicity ( Beers fc Christliebl I2005T ): (4) 
for at least a subset of these CEMP stars, it has been sug- 
gested that their carbon-enhancement reflects the compo- 
sition of th e cloud of gas from which the CE MP star first 
condensed |Rvan et alll2005l ; lAoki et all 120071 ): and (5) the 
observed fraction of CEMP stars in the outer halo compo- 
nent of our Galaxy is roughly twice that of the inner halo 
component, which favours the existence of a high carbon- 
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Figure 10. C and O abundances i n VMP DLAs (red triangles) observed at high spectral resolution. The green triangle represents the 
carbon-enhanced DLA reported by ICooke et al.l ll201ll) . For comparison, we also plot the sample of metal - poor halo stars analysed by 
iFabbian etafl ll2009ah . where we adopt the values for C and O that are calculated assuming the iDrawinl (196SI) formula for collisions 
with hydrogen (Sh = 1; blue squares) We also show a sample of thin- and thick-disc stars (filled and op en circles respectively) from 
IBensbv fc Feltzi ng (2006). The red hatched region corresponds to the transition discriminant outlined by iFrebel. Johnson, fc Bromml 
1120071). where t he un certainty in this relation is shown by the light red shaded region. All of the plotted data have been referred to the 
lAsplund et al.l (2009) solar abundance scale (see Section [SJ. 



producing source ot her than asymptotic giant branch stars 
|Carollo et al.ll201ll ). 

Additional evidence for a high carbon-producing early 
generation of stars has recently been provided by stud- 
ies of the mo st metal-poor DLAs (jPettini et alj 120081; 

iPenprase et al.l |2010T) , which probe entire clouds of near- 
pristine gas. Such surveys find that DLAs and stars tell the 
same story in the metal-poor regime - both show elevated, 
near-solar values of [C/O] when [O/H] ,$ — 2.0. In fact, the 
recent medium spectral resolution study bv lPenprase et all 

(2010) suggests that [C/O] might further increase to super- 
solar values at even lower metallicity, although this still re- 
mains uncertain due to saturation effects. Further evidence 
for an increased carbon yield by an early gener ation of stars 
has re cently come to light with the discovery bv lCooke et al.l 

(2011) of a metal-poor DLA which exhibits a C/Fe ratio 35 
times greater than solar. 

In Figure [TU] we show the updated plot of [C/O] ver- 
sus [O/H] for our full DLA survey, together with values 
for the metal-poor h alo stars (blue squares) analysed by 
IFabbian et al.1 ( 2009a j3, an d a sample of thin- and thick- 
disc stars with C and O abundances determined from 
forbidden lines (filled a nd open black circles respectively; 
IBensbv fc Feltzinell2006h. All data have been corrected for 
the updated lAsplund et al.l (|2009f ) solar abundance scale 
(see Section [5]). We first note that there is generally a 
good agreement - both in the trend and the dispersion of 

9 The data shown in Figure 1101 refer to cross sections f or colli- 
sions w ith hydrogen atoms based on the classical recipe of lDrawm] 
(1969), i.e. using a scaling factor Sh = 1- If hydrogen collisions 
are completely neglected (Sh = 0), this trend is 'stretched' to 
higher [C/O] (by about 0.2 dex) and lower [O/H] (by about 0.3 
dex at the lowest values of [O/H]. 



[C/O] - between the most metal-poor stars and DLAs. The 
new DLA measurements reported here confirm the initial 
indications from th e mo re limited samples co nsidered by 
iPettini et all (|20Qgfl and IPenprase et~ail (|201Ch . The main 
departure fr o m thi s trend is the CEMP DLA reported by 
ICooke et al.l l|201lf) . represented by the green triangle in 
Figure 1101 Aside from this system, no other DLA exhibits 
super-solar [C/O]. This statement is also true for the sam- 
ple of seven 1 ab sorbers a t z abs ~ 6 recently reported by 

iBecker et all (|201ll . see also iBecker et al]|200d). 

It is thus somewhat surprising that IPenprase et al.l 

(|2010h found [C/O] ^ 0.0 for four out of the five VMP DLAs 
in which they could measure this ratio. The difference may 
be due to different sample criteria between our survey and 
theirs: while we have excluded from our analysis absorp- 
tion systems where the Cn lines are saturated, such cases 
may be more difficult to reco gnise at the lower res olution 
of the ESI spectra analysed bv lPenprase et al.l (|201Ch . VMP 
DLAs with super-solar C/O ratios may well exist (and in- 
deed the CEMP DLA discovered by ICooke et all (|201ll ) is 
one such example) , but their C abundance is generally more 
difficult to measure with confidence due to line saturation. 
On the other hand, our survey is not biased against uncov- 
ering systems with lower [C/O] values than those reported 
here; thus, the VMP DLAs in our current sample define a 
lower-envelope in the [C/O] versus [O/H] plane. This en- 
velope appears to be in good agreement with the envelope 
defined by VMP stars (see Figure [TOjl . 

In Figure [TUJ we also show the 'Frebel criter ion' (red 
hatched region; IFrebel. Johnson, fc Bromml 120071 ; see also 
iBromm fc Loebl2003l) which states that, if the fine-structure 
lines of 1 and C 11 dominate the cooling in a near-pristine 
cloud of gas that has been enriched to some critical metal- 
licity, then the first low mass Population II stars will form. 
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Thus, given these conditions, no Population II star should be 
observed in this red hatched region (where the uncertainty 
in this region is given by the light red shaded band). DLAs, 
on the other hand, are not restricted by this criterion. In- 
deed, if a cloud of gas was to be observed in this Population 
II star 'forbidden zone', it may very well form a collection 
of massive stars below or near the critical metallicity! Such 
systems, if found, would provide a unique window to study 
the transition from Population III to Population II star for- 
mation. 

With these considerations in mind, we note that all 
DLAs in our sample will, perhaps unsurprisingly, form Pop- 
ulation II stars. One might, therefore, be tempted to con- 
clude that the stars in the halo of our Milky Way represent 
the same distribution that defines VMP DLAs. To test this 
possibility, we performed a linear fit to [C/O] versus [O/H] 
for the halo stars that have [O/H] ^ —2.0, and calculated 
the deviations about this line for both stars and DLAs. A 
Kolmogorov-Smirnov (K-S) test between the calculated de- 
viations reveals a 75% chance that both VMP stars and 
DLAs are drawn from the same population, which is incon- 
clusive given the present statistics. For this test we have used 
the C/O values in halo stars that were derived assuming ef- 
ficient hydrogen collisions (with a scaling factor Sh = 1; see 
footnote 9). If instead we adopt the stellar values derived 
for inefficient collisions (Sh = 0), we find a 90% chance that 
both samples are drawn from the same population. Such 
good agreement between h alo stars and DLAs d oes no t sup- 
port the recent claim by iTsuiimoto fe Bekkil |201llL who 
suggest that the initial mass function (IMF) of the stars 
that enriched metal-poor DLAs is different from the IMF of 
the stars that enriched Galactic halo stars with their met- 
als. The good general agreement we find between stars and 
DLAs - both exhibiting an elevated C/O ratio at the lowest 
metallicities probed - points to a universal origin for their 
C/O 'excess' in this regime. 



6.2 The O/Fe debate in the metal-poor regime 

We now turn to the relative abundances of oxygen and iron 
at low metallicity. Whilst we cannot do justice to the exten- 
sive literature on this topic, we outline below the basic facts 
that are relevant to our discussion, and direct the interest ed 
reader to the comprehensive review bv lMcWilliaml (| 19971 ). 

The largest oxygen yield comes from the most massive 
stars that explode as SNe II. The budget for iron, on the 
other hand, is largely contribut ed by SNe la w hich typically 
explode ~ 1 Gyr later (see e.g. lGreggicll2010h FI Therefore, 
at early times (when the metallicity is low), one expects O to 
be enhanced relative to Fe. At later times, when the delayed 
contribution of Fe from SNe la kicks in, there is a break in 
the O/Fe trend which is then expected to decrease. Thus, 
the relative abundance of O and Fe allows us to measure 
the relative contribution of SNe la and SNe II (see e.g. th e 
qualitative discussion bv lWheeler. Sneden. fe Truranl fl989). 

In the Milky Way, the break in [O/Fe] occurs roughly at 



10 Iron may also be contributed by SNe la that 'promptly' 
explode at early times (~ 0.1 Gyrs). For the relevant de- 
tails, we direct the interested r eader to the discussion by 
iMannucci. Delia Valle. fc Panagial ^0061 ). 



[Fe/H] ~ —1.0. Whilst there is sound agreement regarding 
the nature of the trend in [O/Fe] when [Fe/H] £ — 1.0, the 
behaviour of [O/Fe] when [Fe/H] ^ — 1.0 is less certain. 
This disagreement stems from the uncertainty of the oxygen 
abundances measured in metal-poor stars: there are four 
different indicators of the oxygen abundance, and to some 
extent t hey all disagree with one another in the metal-poor 
regime (G arci'a Perez et alll2006r i. 

Perhaps the most reliable [O /H] indicator at low metal- 
licity is the forbidden [O i] A6300 line which, despite being 
subject to 3D corre c tions of ~ —0.2 dex when [Fe/H] ~ —2. 
l|Nissen et al.ll2002l ; ICollet. Asplund fe Trampedacbl |2007m , 
is known to form in LTE (|Asplundll2005l ). Unfortunately, 
this line becomes very weak when [Fe/H] < —2.0 and its 
detection requires data of high S/N. After accounting for 
3D corrections to the O abundance, most authors conclude 
that the O/Fe ratio is ap proximately constant at [O/Fe] 
~ +0 .4 for [Fe/H] < - 1 ilNissen et all 120021 ; ICavrel et all 
2004; G arcia Perez et all 120061 ). with perhaps a slight in- 
crease towards the lowest metallicities. 

The most commonly used diagnostic for measuring 
[O/H] in stars is the Ol triplet near 777 nm (A = 7771.9, 
7774.2, 7775.4 A), despite t he fact that i t suffe rs from 
large non- LTE corrections (|Fabbian et all l2009bft . How- 
ever, contrary to the nearly constant value of [O/Fe] below 
[Fe/H] < -1 deduced from the [O i] A6300 line, an LTE anal- 
ysis of the O I triplet leads to a quasi-lin ear increase in [O/Fe] 
with decreasing metallicity (see e.g. iFulbright fe Johnson! 
2003). This discrepancy is often blamed on the uncertain 
(negative) non- LTE corrections to the O I tr iplet. In order 
to mak e headway with the [O/Fe] conflict, iFabbian et all 
(2009b) performed a detailed non- LTE analysis of the O I 
triplet, and found corrections amounting to £ 0.5 dex when 
[Fe/ H] = —3.0, increasing rapidly at lower metallicities (see 
also [Fabbian ct al. 2009 a|). Aft er accounting for such cor- 
rections, Fabbi an et al. (2009b) concluded that almost al 



diagnostics are now conceivably consistent, and that [O/Fe] 
exhibits a roughly fiat plateau with values between +0.4 and 
+0.6 when [Fe/H] < - 1. 

In contrast to the profusion of stellar studies of [O/Fc] 
in the VMP regime, this ratio has received relatively little 
attention in DLAs so far. The reason is that the most readily 
available O I absorption lines are almost always saturated, 
and the weakest lines are of ten blended with unrelate d ab- 
sorption in the Lya forest |Prochaska fe Wolfel I2002I ). For 
these reasons, several authors have investigated the [O/Fe] 
trend in sub-DLAs, where the O I absorptio n lines are weaker 
jPeroux et al.ll2003bl ; IO'Meara et al.ll2005l ). However, uncer- 
tain negative ionization corrections to the Fell lines might 
become important for such systems, complicating the inter- 
pretation. Other authors have instead used [S/Zn] as a proxy 
for [O/Fe] (|Nissen et al.ll2007T l, but these lines become too 
weak in the VMP regime. 

In fact, there have only been two studies in the lit- 
erature th at consider [O/Fe] in DLAs. The first was con- 
ducted by iPetitiean. Ledoux. fe Srianandl (|2008l ). who re- 
ported an [O/Fe] plateau of +0.32 ± 0.10 from their sam- 
ple of 13 DLAs with [Fc/H] < -1.0, (three of which have 
[Fe/H] < —2.0). The second, m ore recent, study was con- 
ducted bv lPenprase et all (|2010T ) whose sample includes five 
DLAs with [Fe/H] < —2.0. Their measurements, however, 
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Figure 11. The [O/Fe] ratio in VMP DLAs (filled red triangle s) whe r e [O/ H] has been measured from high resolution spectra. The green 
triangle refers to the carbon-enhanced DLA reported bv lCooke et al . (2011) ; the upper limit on [Fc/H] and corresponding lower limit on 
[O/Fe] (red triangle with arrows) is for the DLA tow ards J0311— 1722 wh ere Fell absorption is not de tected. The V MP DLA sample is 
compared with the stellar abundance measurements bv lNissen et al] l|2002l ) (circles). ICavrel et al]|2004 ) (squares) and[G arcfa Perez et al] 
(2006) (triangles), all based on the [Ol] A6300 line and corrected for 3D effects (see text for further details). The error bars at the top right 
corner of the plot indicate the typical errors in the stellar abundance measurements. All the measurements in this plot are tabulated, 
together with their individual errors, in Appendix B. 



have large uncertainty (~ ±0.4 dex), so it is difficult to dis- 
cern the underlying trend. 

Our survey constitutes the largest sample of high res- 
olution measures of O I and Fe n absorption in DLAs. The 
[O/Fe] values are plotted in Figure[TT]where, for comparison, 
we also show a selection of [O/Fe] measuremen ts in Galactic 
stars based on the forb i dden [O i] A6300 line jNissen et al] 
120021 ; ICavrel et al1l2004 iGarcia Perez et all l200rj ) with 3D 
corrections applied as we now describe. 

The [0 1] A6300 line corresponds to a forbidden transi- 
tion between two levels of the ground configuration of the 
O I atom, which are closely coupled via collisions. Because 
nearly all oxygen atoms are in the ground state in the at- 
mospheres of late-type stars, one expects LTE to prevail; 
this is confirmed by d etailed statistical equilibrium calcula- 
tions (|Kiselmanll 19931 ). Non-LTE effects on the derived iron 
abundances are also negligible, when lines from the dominat- 
ing io nization stage (Fe n) are considered iMashonkin a et al] 
l201ll ). The 3D - ID corrections are, however, significant 
for metal-poor stars; [O/H] derived from the [Ol] A6300 line 
decreases and [Fe/H] from Fell l ines increases slightly. As 
calculated by iNissen et all l|2002t) . the net effect on [O/Fe] 
for metal-poor main-sequence stars can be approximated by 
the expression [O/Fe] 3D— [O/Fe] id = 0.11 [Fe/H ] , whe reas 
[Fe/H] 3D- [Fe/H] id = -0.04 [Fe/H]. ICavrel et al.l (|2004h as- 
sumed that the same corrections are also valid for metal- 
poor red giants. We have verified that this is approximately 
correct by applying 3D corrections for giant star s as cal- 
culated bv ICollet. Asplund. fc Trampedachl J2007I) for th e 
[0 1] A6300 line and the Fe lines used bv lCavrel et al l J2004h . 
The same 3D corrections are then also expected for the 
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Figure 12. The distribution of [O/Fe] values (where [Fe/H] ^ 
—2.0) in DLAs for our sample (red histogram) as compared to the 
compilation of [O/Fe] in stars measured from the [Ol] A6300 line 
(blue histogram). The dashed vertical red and blue lines indicate 
the median values for DLAs and stars respectively, corresponding 
to a difference of only ~ 0.01 dex. 

cool subgiants studied by iGarcfa Perez et al] ()2006l L be- 
cause they have atmospheric parameters intermediate be- 
tween those of the main-sequence and red giant stars. For 
reference, we list in Appendix[B]values of [O/Fe] and [Fe/H] 
corrected for 3D effects, as well as the [O/Fe] and [Fe/H] 
values measured in our sample of VMP DLAs. 

Once the above-mentioned 3D corrections are applied, 
it can be seen that stars and DLAs share a similar trend of 
[O/Fe] with decreasing metallicity. To illustrate this point, 
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Table 13. The mean value and dispersion in X/O for each 

ELEMENT OF THE VMP DLAS IN OUR SAMPLE. 



X 


[<x/o>] 


CT [(x/o>] 


«[(X/o>] 


c 


-0.28 


0.12 


9 


N 


-1.05 


0.19 


13 


Al 


-0.44 


0.13 


11 


Si 


-0.08 


0.10 


21 


Fe 


-0.39 


0.12 


20 



we show the two samples (where [Fe/H] ^ —2.0) as his- 
tograms in Figure 1121 A K-S tes|^j reveals that the prob- 
ability the two data sets are drawn from the same parent 
population is 71%. 

The DLA values of [O/Fe] exhibit relatively little scat- 
ter given the errors. In the range —3 ^ [Fe/H] —2, [O/Fe] 
is consistent with a constant value: [(O/Fe)] = +0.35 ±0.09. 
This is in good agreement with the m ean value reported by 
iPetitiean. Ledoux, fe Srianandl |200Sl ). [(O/Fe)] = +0.32 ± 
0.10, for -2.0 £ [Fe/H] < - 1.0. Interestingly, there may 
be a hint in our data that [O/Fe] increases further when 
[Fe/H] ;$ — 3 (see Figure [TTj) . which would presumably be 
indicative of a contribution from more massive stars. How- 
ever, more data are required to confirm this 'trend' which 
at present is suggested by the two most metal-poor DLAs in 
our sample. For the moment, we simply conclude that the 
'cosmic' trend of [O/Fe] in the VMP regime (-3 ^ [Fe/H] 
^ —2) reaches a plateau of ~ +0.35, and is remarkably tight, 
especially given the observational errors. 

In closing, the DLA measurements of [O/Fe] help re- 
solve the controversy regarding the relative abundances of O 
and Fe metal-poor Galactic stars. It is plausible that VMP 
DLAs harbour the reservoir of neutral gas that will later 
condense to form a population of VMP stars. It is thus ex- 
pected that towards the lowest metallicities both stars and 
DLAs should exhibit a similar trend. We conclude that, un- 
less there are marked differences between the chemical evolu- 
tion his tories of DLAs and the early Ga l axy, o ur results and 
those of lPetitiean. Ledoux. fe Srianandl C 2008) favour an ap- 
proximately constant plateau of stellar [O/Fe] values when 
[Fe/H] ;$ — 1, with perhaps a mild increase with decreasing 
[Fe/H]. In any case, given the on-going improvement in the 
accuracy of the stellar models and the increasing samples 
of VMP DLAs, we anticipate that this issue may well be 
settled in the near future. 



7 DISCUSSION 

7.1 The typical VMP DLA 

With the large sample of measurements we have assembled, 
we can now attempt to reconstruct the abundance pattern 
of a 'typical' VMP DLA, and to consider the clues it may 
provide on the nucleosynthesis by the earliest generation of 
stars. To achieve this goal, we are obviously required to se- 
lect a reference element other that H, as we have no means to 

11 We have not in cluded the star CS 22949-037 in this test (from 
ICavrel et al ■I2004IV since it exhibits a peculiar abundance pattern, 
with [O/Fe] = +1.54. 
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Figure 13. The abundance pattern of a typical VMP DLA is 
illustrated by the black boxes, where the height in each box rep- 
resents the dispersion in the population. The dashed line corre- 
sponds to the solar abundance ratios. 



determine how much H was mixed with the nucleosynthetic 
products from the earliest generations of stars. Rather, the 
ratio of two metals provides the best handle for determining 
the properties of the generation of stars from which they 
were synthesised. 

With our goal to probe early nucleosynthesis borne in 
mind, the most appropriate reference element is O, for the 
following reasons: (1) the dominant O yield comes from a 
single source - massive stars. Thus, the origin of O is well- 
understood; (2) O is the most abundant metal in the Uni- 
verse; and (3) at the lowest metallicities, where we expect 
to uncover the signature from early nucleosynthesis, several 
O lines become unsaturated in DLAs. Thus, it is relatively 
straightforward to measure [O/H]. Taking oxygen then as 
the reference element, we have constructed the typical abun- 
dance pattern for a VMP DLA by determining the mean 
(X/O) ratio for each available element, X, and then referring 
this mean value to the adopted solar scale (i.e. we take the 
log of the mean, [(X/O)]). These mean values are listed in 
Table [131 along with the dispersion in the available measure- 
ments (cr^x/o)])- In the last column of Table [1U1 we indicate 
the total number of DLAs that were used in determining the 
mean X/O. 

The corresponding abundance pattern is illustrated in 
Figure 1131 In the following subsection, we investigate the 
most likely origin of the metals in VMP DLAs by directly 
comparing this typical abundance pattern to model yield 
calculations of both Population II and Population III stars. 
Before continuing, however, it is important to keep in mind 
that the 'typical' [(N/O)] in Table HUl may be biased high, 
because it does not include a number of upper limits, where 
the N I lines are too weak to be measured. As discussed above 
(Section [STTJ, it is also possible that [(C/O)] may be biased 
towards lower values than the true mean. 

According to our sample definition (see Section I4.8f) . 
this 'typical' abundance pattern is based on DLAs where 
[Fe/H] ^ —2.0. By restricting this sample to only contain 
DLAs with [Fe/H] ^ —2.5, we found that the typical abun- 
dance for all elements considered here changes by no more 
than 0.02 dex, except in the case of [(N/O)], which is lower 
by 0.14 dex. Thus, our choice of metallicity cut does not in- 
troduce a significant bias into the typical VMP DLA abun- 
dance pattern. 
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7.2 Clues to early episodes of nucleosynthesis 

There is increasing evidence to suggest that the most metal- 
poor DLAs may retain th e signature from the earliest 
episodes of star format ion |Erni et al.l 120061 ; iPettini et al.l 
120081 ; ICooke et al.ll201lh . In this picture, the most metal- 
poor DLAs condensed directly out of material that was 
enriched by either: (1) an external halo that distributed 
its pro ducts over large cosmological v olumes via multi-SN 
events l)Madau. Ferrara. fc Reeal200ll L or perhaps (2) just 
a few SNe from the halo in which the DLA now resides 
(|Bland-Hawthorn et afll201ll ). 

Cosmological simulations of galaxy formation support 
the possibility that such DLAs sti ll retain the chem- 
ical signature of early enrichment (|Pontzen et al.l 120081; 
iTescari et aT1 [2009); the most metal-poor DLAs arise in low 
mass halos that have undergone little to no in situ star for- 
mation. It is possible, however, that VMP DLAs acquired 
some of their metals at later times from nearby sources that 
deli vered metals into the IGM via galactic superwinds (see 
e.g. lOppenheimer fc Dave1l2008h , which may complicate the 
interpretation. Perhaps the most straightforward way to dis- 
criminate between these enrichment scenarios is to compare 
the model yields of both Population III and Population II 
stars with that of the typical VMP DLA described in Sec- 
tion EH 

We consider three sources that could be responsible 
for the metals in VMP DLAs: (1) massive metal-free stars, 
with main sequence masses in the range 1 40 — 260 Mq that 
explo de as pair-instability SNe (PISN; Irleger fc Wooslevl 
|2002| ) ; (2) massive metal-free stars, with progenitor masses 
in the r ange 10 — 100 Mq that e xplode as core-collapse SNe 
fCCSN: iHeger fc Wooslevl l201Ch ; and (3) massive Popula- 
tion II (and I) stars, with progenitor mass es in the range 13— 
35 M e, covering a range in metallicity l)Chieffi fc Limongil 
2004), also ending their lives as CCSN. To determine the 
dominant source of the metals in VMP DLAs, we inte- 
grate these model yields over a Salpeter-like power law IMF, 
dN/dM <x M" (1+7) (where 7 = 1.35 for a Salpeter IMF), 
and consider three values for the power law index in the 
case of zero metallicity (7 = 1.35,2.35,3.35). For massive 
Population II stars, we consider only a Salpeter IMF, with 
7 = 1.35. The results of these calculations are shown in 
Figure Q3] 

Let us first consider the yields from metal-free stars 
with masses in the range 140 — 260 Mq. These stars ex- 
plode as PISN, the physics of which is well-u nderstood 
jHeger fc Wooslevll2002l : lUmeda fc Nomotdl2002l ). The cal- 
culated yields from PISN are thus the least model dependent 
of the three cases considered here. Qualitatively, these mod- 
els could have been ruled out on the basis of the near-solar 
[Si/O] that is typical of the VMP DLA population; PISN are 
expected to yield supersolar [Si, S, Ar, Ca/O]. Indeed, for 
the range of 7 considered here, such SNe provide a poor fit to 
the typical VMP DLA population (top panel of Figure [14"]) . 

We now turn to models of massive stars (10 — 100 Mq) 
that end their lives as CCSN. Unfortunately, the explosion 
mechanism of CCSN is poorly understood, and several un- 
known physical effects need to be parameterized and suit- 
ably adjusted to find the best solution for a given set of 
data. In particular, one usually parameterizes the explosion 
energy, the degree of mixing between the stellar layers dur- 
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Figure 14. The abundance pattern of the 'typical' VMP DLA 
(black boxes; cf. Figure 113)1 is compared to nucleosynthesis mod- 
els of massive stars (symbols connected by lines). Top panel: 
The symbols represent the yield s from pair-instability su pernovae 
(PISN) of zero metallicity stars (Hcgcr & Woosley 2002) for three 
indices of a power law IMF (red, green, blue corresponds to 
7 = 1.35,2.35,3.35 respectively). Middle panel: Same as above, 
except for c ore-collapse supernovae ( CCSN) models of zero metal- 
licity stars l|Heger fc Wooslevll20icf ). Bottom panel: Comparing 
the typical DLA abundance patter n with the model explo sive 
yields of massive Population II stars l|Chieffi fc Limongill2004h . A 
range of metallicities is considered, as indicated by the accom- 
panying legend. In all cases for this bottom panel, we plot the 
integrated yields from the Salpeter IMF with a power law index 
of 7 = 1.35. 



ing the explosion, and the amount of material that falls back 
onto the central remnant. 

The most recent suite of published CCSN yi elds for 
metal-free stars are those bv lHeger fc Wooslevl l|2010h . These 
computations provide a detailed account of the nucleosyn- 
thetic products over a large range of progenitor masses 
(10 — 100 Mq) with a typical mass resolution of 0.1 Mq. 
As part of their study, Irleger fc Wooslevl (|2010n compiled a 
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database of these model yields that are imported into their 
STARFllQ software. This software is designed to objectively 
sieve through the vast parameter space and select the explo- 
sion parameters that best fit the data. 

To maintain consistency with the other yield models 
that are considered here, we compare the abundance pat- 
tern of the typical VMP DLA to the expected yields for 
the three power law indices of a Salpeter-like IMF (7 = 
1.35,2.35,3.35). We therefore froze the remaining parame- 
ters to the 'standard' case, which corresponds to a constant 
explosion energy (1.2 x 10 51 erg) for all masses in the range 
10-100 M (see lHeger fc WooslevllioTol for further details). 
This standard case sets the piston location of the explosion 
to be at the base of the oxygen burning shell (where the 
entropy per baryon ~ 4), and applies a mixing boxcar filter 
with a width which is 10% the He core size. The material 
that falls back onto the central remnant is not parameter- 
ized in this code, but is instead calculated by switching the 
piston off 100 s after the explosion and defining an inner 
boundary condition where material is accreted. 

All three fits are reproduced in the middle panel of Fig- 
ure 03] The 'standard' case with the yields from metal- free 
stars seems to produce a reasonable agreement with the ob- 
served metal ratios in the 'typical' VMP DLAs, although 
Al is discrepant by ~ 0.5 dex and, interestingly, IMF slopes 
steeper than Salpeter seem to fit the data best. Of course, by 
relaxing some of the default constraints it may be possible 
to improve the fit further, but we have refrained from doing 
so, given that there are still many uncertainties in accurately 
modelling the physics behind the explosion. Some of these 
uncertainties are only now beginning to be addresse d in 
some detail (see iJoggerst. Almgren. fc Wooslevl l2010bl and 
references therein). 

Finall y, we consider the set of model CCSN yields pub- 
lished by IChieffi fc Limongil <|2004h . which allow us to test 
whether or not Population II stars can also account for the 
origin of the metals in VMP DLAs. Before comparing the 
models by Chieffi & Limongi to the typical VMP DLA, it 
is worth noting the important diff erences between this cod e 
and the one described above by iHeger fc Wooslevl (|2010h . 
As ide from the obvious diffe rence in metallicity, the models 
bv lChieffi fc Limongil (|2004h target the mass range 13 — 35 
Mq with a relatively coarser mass resolution of 5 Mq. In ad- 
dition, this code parameterizes the amount of material that 
falls back onto the central remnant. In their standard case, 
this prescription requires 0.1 Mq of 56 Ni to be ejected from 
the star, and thus all material interior to this mass coordi- 
nate is 'accreted' by the remnant. Finally, this code is yet to 
implement a 'mixing parameter' to model the mixing that 
takes place between the stellar layers during the SN explo- 
sion. Adopting their standard case, which has an explosion 
energy of 1.2 x 10 51 erg, we integrate the model yields over 
a Salpeter-like IMF with 7 = 1.35. The results are shown in 
the bottom panel of Figure 1141 

These calculations adopt an initial chemical composi- 
tion that is s imply scaled from the sola r abundance pattern. 
According to lChieffi fc Limongil (j2004h . the model yields for 



12 STARFIT is written with the INTERACTIVE DATA LANGUAGE soft- 
ware and is available from: 
http://homepagcs.spa.umn.edu/~alex/znuc/ 



stars with an initial metallicity of Z ^ 0.005 Zq are not 
strongly dependent on the initial composition of the star. 
It is important to note, however, that by introducing an 
a-enhancement to the initial metallicity (which is perhaps 
more realistic than simply scaling the solar abundance pat- 
tern), the yields for the odd atomic number elements are 
increased (see their Figure 1). Thus, our calculations may 
underestimate the N/O and Al/O ratios. Furthermore, at 
higher metallicities (Z > 0.005 Zq), the yields do depend on 
the initial composition of the stars. 

Inspecting the bottom panel of Figure 1141 it can been 
seen that, at face value, the Z/Z© = 0.0, 0.05 and 0.3 mod- 
els provide reasonable fits to the abundance pattern of a 
typical VMP DLA. We also note the broad agreem e nt be - 
twe en the metal-free model s by IHeger fc Wooslevl (|2010h 
and lChieffi fc Limongil \2004 ). On the other hand, the Z/Z 
= 0.05 and 0.3 models also seem to provide reasonable fits 
to the abundance pattern. However, once the stars have 
reached metallicities greater than 1/20 of solar, the metal 
ratios in the gas should be compared with the predictions of 
full chemical evolution models which are beyond the scope 
of this paper. 

Given the current (largely) model-dependent nature of 
these calculations, we are unable to draw firm conclusions at 
this stage. Whilst the above models suggest that metal-free 
stars could have synthesised the metals that now reside in 
VMP DLAs, we cannot rule out the possibility that Popu- 
lation II stars are also responsible. We suspect that it will 
be necessary to measure the abundances of additional met- 
als in order to better distinguish between Population II and 
Population III models, since the ratios of the most abun- 
dant elements ([C/O], [Si/O] and [Fe/O]) provide the weak- 
est constraints on the nature of the objects that synthesised 
them. Some of the largest differences between these models 
are exhibited by the iron-peak elements, and in particular, 
Ti, Ni and Zn. The detection of these rarer elements in the 
most metal-poor DLAs will have to wait for the advent of 
the next generation of extremely large optical telescopes. 

Finally, we note that the models used here to com- 
pare with a typical VMP DLA are still quite dependent 
on unknown physics; the single largest uncertainty in these 
models is the explosion mechanism. Additional physics also 
needs to be included, such as the mixing induced by stellar 
rotation |Meynet. Ekstrom. fc Maederl l200e3 : iHirschil I2007I : 
iMevnet et al.ll2~010l;|joggerst et al]l2010af) and the R avleigh- 
Taylor instability (jJoggerst. Wooslev. fc Hegerll2009l ). These 
effects will presumably be considered in the next generation 
of fine-grid nucleosynthesis models, when the limitations of 
computing power will hopefully be less of a concern. 



7.3 Comparison with data of medium spectral 
resolution 

Finally, we compare DLA abundance determinations ob- 
tained from high (R ~ 40000) and medium (R ~ 5000) 
spectral resolution data. Such a comparison is motivated by 
the realization that, even with efficient echelle spectrographs 
on 8-10 m telescopes, it is typically necessary to integrate 
on a single QSO for the equivalent of one night in order to 
obtain the S/N ratio required to measure elemental abun- 
dances from high resolution spectra. By settling for lower 
resolutions, the exposure times are greatly reduced; for ex- 
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Table 14. COLUMN DENSITIES ESTIMATED FROM HIGH AND MEDIUM 
SPECTRAL RESOLUTION SPECTRA 



Ion 


log 7V a log N b 

6 J "coo o pcn.m 


log N£ enc 






J0831+3358 




Al II 


12 


.19 ±0.06 


> 11.56 


Sin 


13. 


.75 ± 0.04 


> 13.03 


Fe II 


13. 


.33 ± 0.06 


> 13.13 






J1001+0343 




Cn 


13. 


.58 ±0.02 13.63 ±0.06 


13.76 ±0.13 


Oi 


14. 


.25 ±0.02 13.90 ±0.07 


13.98 ±0.08 


Sin 


12 


.86 ±0.01 12.70 ±0.05 


12.82 ±0.12 






J1037+0139 




Oi 


15 


.06 ± 0.04 


> 14.48 


Alii 


12 


.32 ±0.03 


> 12.09 


Sin 


13. 


.97 ±0.03 


> 13.69 


Fen 


13. 


.53 ±0.02 


> 13.33 



a Our results. 

b Column densities measured bvlPenp rase et aL I d2010f) . 
c Column density estimates bv lPenprase et aL I" il2010h after 
applying a saturation correction. 

ample , most of the QSOs in the survey by iPenprase et al.l 
(2010) were observed for about one hour with the Echelle 
Spectrograph and Imager (ESI) . Clearly, it is of interest to 
test how similar the abundance measurements are for the 
most metal-poor DLAs, if we were to forgo the accuracy of 
high spectral resolution in order to secure a larger sample. 

In this context, there are two main concerns that po- 
tentially limit the accuracy of abundance measurements 
from medium (as opposed to high) spectral resolution data. 
First, VMP DLAs t y pically have line wid t hs ^ 10 km s" 1 
(iLedoux et all l200rj ; iMurphv et""aH 120071 ; iProchaska et all 
2008), which are unresolved at R — 5000. One must there- 
fore appeal to a curve-of-growth analysis appropriate to a 
single absorbing cloud, w hich is often an oversimplification 
(see our Table [1 and also IProcha ska 2006). Second, the rel- 
evant absorption lines may be saturated, and the degree of 
saturation may be difficult to estimate correctly, even with 
a well-defined curve-of-growth. 

As it happens, three VMP DLAs from the present 
work are in common with the lower resolution survey of 
IPenprase et alj (|2010l ). providing us with the means to com- 
pare column densities from the two sets of spectra, as in Ta- 
blgJTg] For two of the DLAs, J0831+3358 and J1037+0139, 
IPenprase et alj |2010l ) reported lower limits on the column 
densities of the available metal ions. As can be seen from Ta- 
ble !14l while these lower limits are always consistent with the 
values measured from our echelle spectra, they fall short of 
the true column density by widely differing amounts, from as 
little as —0.13 dex to as much as —0.72 dex. This wide range 
significantly reduces the usefulness of the lower limits. 

Turning now to the DLA at 2 a b s = 3.07841 towards 
J1001+0343, we recall that our UVES spectra show that 
the absorption arises in a single component with Doppler 
parameter b = 7.0 ± 0.1 km s -1 (see Figure [3|. This value 
is not too diss i milar from b — 7.5 km s _1 estimated by 
IPenprase et"afl (|2010l ) from a curve-of-growth analysis. In- 
deed, the two analyses give consistent estimates of the Sin 



column density (after [Penprase et al.l l|2010l ) apply a satura- 
tion correction). The Cn a nd O I lines, however, tell a dif- 
ferent story. For this DLA, IPenprase et alj (|2010T ) need not 
have applied a saturation correction to the CllA1334 line, 
since it is not strongly saturated (see Figure(3|. Indeed, prior 
to applying such a correction, their column density estimate 
was in broad agreement with that derived here. Conversely, 
OlA1302 is clos er to saturation and, even with the correc- 
tion applied by IPenprase et all (|20ld ). these authors' esti- 
mate falls short of the value deduced here by nearly a factor 
of two. The combined effect is an overestimate of [C/O] by 
0.45 dex. 

Based on this example, there appear to be non- 
negligible uncertainties in the derivation of element ratios 
from spectra at R ~ 5000. While these u ncertainties did 
not af fect the principle goal of the study bv lPenprase et all 
(2010) to uncover the most metal-poor DLAs - it would 
appear that high resolution observations are indeed neces- 
sary to measure element abundances in VMP DLAs with an 
accuracy better than a factor of ~ 2. 



8 SUMMARY AND CONCLUSIONS 

We have conducted a survey for very metal-poor DLAs to 
shed light on the earliest episodes of nucleosynthesis in our 
Universe. Our sample includes seven new DLAs observed 
with high resolution spectrographs (R Si 30000); when com- 
bined with the five met al-poor DLAs previously reported 
from this programme jPettini et al.l 120081 ; I Cooke et al.1 
l201ll ) and an additional ten DLAs from the literature, it 
constitutes the largest survey to date for DLAs with a 
metallicity [Fe/H] < —2.0. From the analysis of these data, 
we draw the following conclusions. 

(i) Having now doubled the sample of DLAs where the 
C/O ratio is measured from unsaturated absorption lines, 
we confirm that DLAs exhibit near-solar values of C/O at 
the lowest metallicities probed. Furthermore, we find good 
agreement in the C/O ratio observed in our sample of DLAs 
and in recent compilations of the most metal-poor Galactic 
halo stars. We argue that such good agreement points to a 
universal origin for the C/O 'excess' in this regime. 

(ii) For the first time, we investigate the [O/Fe] ratio in 
very metal-poor DLAs. For 20 DLAs with [Fe/H] < -2.0, 
we find a small dispersion around a mean value [(O/Fe)] = 
+0.39 ±0.12. We have also presented tentative evidence for 
a rise in the [O/Fe] ratio when [Fe/H] ,$ — 3.0. 

(iii) In view of the long-standing debate as to the behaviour 
of the [O/Fe] ratio in metal-poor Galactic halo stars, we have 
compared the stellar trend to that observed in our sample 
of DLAs. We find good agreement between stars and DLAs 
when the stellar oxygen abundance is measured from the 
[Ol]A6300 line (after correcting for 3D effects). Based on 
the available DLA samples, we conclude that [O/Fe] is essen- 
tially flat in the metallicity interval —3.0 [Fe/H] ;$ —1.0, 
with the possibility of an increase at yet lower metallicities. 

(iv) We have constructed the abundance pattern of a typ- 
ical very metal-poor DLA for the five most commonly ob- 
served metals, using O as a reference. We find that Si/O is 
just below solar ([(Si/O)] = -0.08), whilst [(C/O)] = -0.28 



© 2011 RAS, MNRAS OOO.rTJMl 



24 Cooke et al. 



and [(Fe/O)] = —0.39. The largest deviations from a solar 
scaled abundance pattern are exhibited by N and Al, with 
[(N.Al/O)] = -1.05,-0.44. 

(v) One of the main aims of this work was to investigate 
the origin of the metals in the most metal-poor DLAs. To 
achieve this goal, we compared the abundance pattern of a 
'typical' VMP DLA with those expected from model calcu- 
lations using the yields of Population II and Population III 
stars. For the few elements considered here, we find a reason- 
able agreement between the abundance pattern of the typi- 
cal VMP DLA and the 'standard model' of a population of 
metal-free stars (i.e. a top-heavy initial mass function where 
all stars explode as core-collapse supernovae with an energy 
of 1.2 x 10 51 erg). However, given that we only have access 
to a handful of metals, we cannot unambiguously rule out 
(an additional contribution from) more metal-rich Popula- 
tion II stars. On the other hand, we are able to firmly con- 
clude that the typical very metal-poor DLA was not solely 
enriched by pair-instability supernovae from very massive 
metal-free stars. 

Our ongoing programme to measure the abundances 
in the most metal-poor DLAs complements local studies of 
Galactic metal-poor halo stars. The good agreement we have 
found between these two populations suggests a universal 
origin for their metals. The results presented here empha- 
size the importance of measuring elemental abundances in 
the most metal-poor DLAs; these systems present us with a 
unique window of opportunity to probe the nucleosynthesis 
by some of the earliest structures in the Universe. 
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APPENDIX A: VMP DLA COLUMN 
DENSITIES 

To facilitate comparison with future data sets, Table IA1I 
lists ion column densities for the DLAs considered in this 
study. The corresponding element abundances are collected 
in Table [Til 
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Table Al. C, N, O, Al, Si, and Fe column densities in VMP DLAs 



QSO 


^abs 


logAT(Hl) 


logAf(Cn) 


logTV(Ni) 


logTV(Oi) 


logJV(Aln) 


log7V(Si n) 


logA^Fell) 


Ref. a 






(cm -2 ) 


(cm -2 ) 


(cm -2 ) 


(cm -2 ) 


(cm -2 ) 


(cm -2 ) 


(cm -2 ) 




Our Metal-Poor DLA Sample 


















J0035-0918 


2.34010 


20.55 ±0.10 


10.4/ ± U.1D 


iq ri _|_ n flft 

lo.Ol ± U.Uu 


i a nft _i_ n no 
14. yt) ± U.Uo 


1 1 7Q j_ n nc; 
11. /o ± U.UO 


13.41 ±0 


01 


lz.yo ± U.U ( 


z 


J0311-1722 


3.73400 


20.30 ± 0.06 


14. uz it U.Uo 




i /i 7H -U n nc 
14. / U jr. U.Uo 




13.31 ± 


07 


^, lo. (D 


1 


J0831+3358 


2.30364 


20.25 ± 0.15 




^ 12.78 


14.93 ±0.05 


12.19 ±0.06 


13.75 ± 


01 


13.33 ± 0.06 


1,4 


Q0913+072 


2.61843 


20.34 ± 0.04 


13.98 ±0.05 


12.29 ±0.12 


14.63 ±0.01 


11.78 ±0.03 


i q on _|_ n 


01 


12.99 ± 0.01 


3 


J1001+0343 


3.07841 


20.21 ± 0.05 


lo.oo ± U.U2 


^ 12.50 


i/i occ _i_ n no 
14. ZO ± U.UZ 




i 2 86 4- 

J. _ . OU _1_ 


01 


i o c;n _l_ n 1 /t 
12. 5U ± U.14 


1 


J1016+4040 


2.81633 


19.90 ± 0.11 


1 Q ftft _l_ C\ (\A 

lo.ou ± U.U4 


^ 12.76 


i/i iq _i_ n no 
14. lo ± U.Uo 




1 2 qn 4- n 

J. _ . Z)\J _1_ VI 


05 




3 


J1037+0139 


2.70487 


20.50 ± 0.08 




i q orr _L n c\a 
16.2, 1 ± U.U4 


i c: nft _L n n/i 
lo. Uu ± U.U4 


i o q o _i_ n nQ 
IZ.oZ ± U.Uo 


13.97 ± 


03 


1 q c; q _L n no 
lo.OO ± U.U2 


1 


J1340+1106 


2.50792 


20.09 ± 0.05 




i o on _L n c\a 
12. oU ± U.U4 


i c: no _L n no 
15.U2 ± U.Uo 


i o 17 _i_ n no 
1Z.Z / ± U.UZ 


13.75 ± 


02 


1 q aci _i_ n no 
lo.4y ± U.U2 


1 


J1340+1106 


2.79583 


21.00 ± 0.06 




i a c\a _i_ n no 
14. U4 ± U.UZ 


i ft c\a _i_ n n/i 
lu.U4 ± U.U4 


i q o a _i_ n nQ 
lo.24 ± U.Uo 


14.68 ± 


02 


1 a qo _L n m 
14. oz ± U.U1 


1 


J1419+0829 


3.04973 


20.40 ± 0.03 




i q oq _L n no 
lo.zo ± U.UZ 


i c 17 i n no 
10.1 / ± u.uz 




1 3 83 + 


01 


1 q K/t _L n nQ 
lo.04 ± U.Uo 


i 

l 


J1558+4053 


2.55332 


20.30 ± 0.04 


i/i n _l n nft 
LQ.ZZ ± U.Uu 


i o ftft _i_ n Q7 
12. DO ± U.U/ 


1/1 K /I _L C\ (\A 

14.04 ± U.U4 


11 no _i_ n nft 
11. yz ± U.Uu 


1 3 32 4~ 


02 


1 q c\i _i_ n nft 
lo.U/ ± U.Uo 


3 


Q2206-199 


2.07624 


20.43 ±0.04 


14.41 ± U.Uo 


i o 7o i n n^ 
1Z. / y ± U.Uo 


i c; ne: j_ n nQ 
lo.Uo ± U.Uo 


i o ic j. n m 
IZ.lo ± U.U1 


1 S 65 4- 

-LO. UtJ _1_ VJ 


01 


1 q qq _l n ni 
lo.oo ± U.U1 


Q 

o 


Literature 


DLAs 




















Q0000-2620 


3.39012 


21.41 ±0.08 




i a ic\ _i_ n no 
14. / U ± U.Uz 


is /io_Ln in 
1U.4Z ± U.1U 




15.06 ±0 


02 


1 a Q7 _l n nQ 
14. o/ ± U.Uo 





Q0112-306 


2.41844 


20.50 ±0.08 




13.16 ±0.04 


14.95 ±0.08 




13.62 ±0 


02 


13.33 ±0.05 


6 


J0140-0839 


3.69660 


20.75 ±0.15 


14.13 ±0.08 


sC 12.38 


14.69 ±0.01 


11.82 ±0.04 


13.51 ±0 


09 


12.77 ±0.19 b 


7 


J0307-4945 


4.46658 


20.67 ±0.09 




13.57 ±0.12 


15.91 ±0.17 


13.36 ±0.06 


14.68 ±0 


07 


14.21 ± 0.17 


8 


Q1108-077 


3.60767 


20.37 ±0.07 




sC 12.84 


15.37 ±0.03 




14.34 ±0 


02 


13.88 ±0.02 


6 


J1337+3153 


3.16768 


20.41 ±0.15 


13.98 ±0.06 


sC 12.80 


14.43 ± 0.09 


12.00 ±0.05 


13.24 ±0 


05 


13.14 ±0.26 


9 


J1558-0031 


2.70262 


20.67 ±0.05 




14.46 c 


15.86 c 




14.24 c 




14.11 c 


10 


Q1946+7658 


2.84430 


20.27 ±0.06 




12.59 ±0.04 


14.82 ±0.01 




13.60 ±0 


01 


13.24 ±0.01 


11 


Q2059-360 


3.08293 


20.98 ± 0.08 




13.95 ±0.02 


16.09 ± 0.04 




14.86 ±0 


05 


14.48 ± 0.02 


6 


J2155+1358 


4.21244 


19.61 ±0.10 


13.95 ± 0.06 




14.50 ± 0.05 


11.92 ±0.17 


13.25 ± 


04 


12.93 ± 0.23 


12 



References— 1: This work; 2 : [Cooke et all ll201lll; 3: IPettini et al] teOOSl): 4: IPenprase et a I J2010h : 5: iMolaro et al.M200(t) : 6: 
Petitiean. Ledoux, fc Sriaiiandl Impost); 7: lEllison et al] ||201(]|'); 8: iDessauges-Zayadskv et al.l j200j) ; 9: ISrianand et alj ^OlOp ; 10: 
O'Meara et al l j2006l 'l : 11: IProchaska et alj l|2002l 'l; 12: iDessauges-Zavadskv et al.1 (|2003^ 



lEllison et al. fl l2010h quote a 3<r upper limit to the Fell column density of log JV(Fe 11) /cm 2 < 12.73. We have since rereduced these 
data (as described in Section [2}, and detected the FellA1608 line at the 4<r level. The column density for Fell quoted here is derived 
using the optically thin limit approximation. 

c An error estimate for this measurement was not provided by the authors. 



APPENDIX B: OXYGEN AND IRON IN DLAS 
AND STARS 

In Tables [BT] and [H we list the values of [O/Fe] and [Fe/H] 
plotted in Figure QT] for DLAs and Galactic stars respec- 
tively, together with the corresponding errors. The stellar 
measurements are all from the [O i] A6 3 00 line (taken from 
Nissen et ail |2002| ; ICavrel et al] 120041 : iGarcfa Perez et al.1 
2006J), with appropria te 3 D corrections based on the work 
bvlNissen et al) (|2002l ) and lCollet, Asplund. fc Trampedacbl 
(|2007l ); see text in Section f6. 21 for further details. 
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Table Bl. [Fe/H] and [O/Fe] in VMP DLAs 



QSO Name 




^abs 


lc 


>giV(Hi) 
(cm" 2 ) 


[Fe/H] 


[O/Fe] 


Ref. a 


Q0000— 2620 


3. 


.39012 


21 


.41 ± 0.08 


—2.01 ± 0.09 


+0.33 ± 0.10 


2 


J0035— 0918 


2 


.34010 


20 


.55 ± 0.10 


—3.04 ± 0.12 


+0.76 ± 0.11 


3 


Q0112— 306 


2 


.41844 


20 


.50 ± 0.08 


—2.64 ± 0.09 


+0.40 ± 0.09 


1 


TA I Aft flOOA 

J0140— 0839 


3 


.69660 


20 


.75 ± 0.15 


—3.45 ± 0.24 


+0.70 ± 0.19 


5 


J 0307— 4945 


4 


.46658 


20 


.67 ± 0.09 


— 1.93 ± 0.19 


+0.48 ± 0.24 


6 


J0311 — 1722 


3. 


.73400 


20 


.30 ± 0.06 


^ —2.03 


> —0.26 


1 


J0831+3358 


2 


.30364 


20 


.25 ± 0.15 


—2.39 ± 0.16 


+0.38 ± 0.08 


1.7 


Q0913+072 


2 


.61843 


20 


.34 ± 0.04 


—2.82 ± 0.04 


+0.42 ± 0.02 


8 


J1001+0343 


3 


.07841 


20 


.21 ± 0.05 


—3.18 ± 0.15 


+0.53 ± 0.14 


1 


T "i noT i *a "i or\ 

J 1037+0139 


2 


.70487 


20 


.50 ± 0.08 


r\ A a 1 r\ c\a 

—2. 44 ± 0.08 


+0.31 ± 0.04 


1 


Q1108— 077 


3. 


.60767 


20 


.37 ± 0.07 


— 1.96 ± 0.07 


+0.27 ± 0.04 


i 


J1337+3153 


3 


.16768 


20. 


.41 ±0.15 


-2.74 ±0.30 


+0.07 + 0.28 


9 


J1340+1106 


2 


.50792 


20. 


.09 ±0.05 


-2.07 ±0.05 


+0.31 + 0.04 


1 


J1340+1106 


2 


.79583 


21. 


.00 ± 0.06 


-2.15 ±0.06 


+0.50 + 0.04 


1 


J1419+0829 


3 


.04973 


20. 


.40 ± 0.03 


-2.33 ±0.04 


+0.41 + 0.04 


1 


J1558-0031 


2 


.70262 


20. 


.67 ±0.05 


-2.03 b 


+0.53 b 


10 


J1558+4053 


2 


.55332 


20. 


.30 ±0.04 


-2.70 ±0.07 


+0.25 + 0.07 


8 


Q1946+7658 


2 


.84430 


20. 


.27 ±0.06 


-2.50 ±0.06 


+0.36 + 0.02 


11 


Q2059-360 


3 


.08293 


20. 


.98 ±0.08 


-1.97 ±0.08 


+0.39 + 0.04 


4 


J2155+1358 


4. 


.21244 


19. 


.61 ±0.10 


-2.15 ±0.25 


+0.35 + 0.24 


12 


Q2206-199 


2 


.07624 


20. 


.43 ± 0.04 


-2.57 ±0.04 


+0.50 + 0.03 


8 



^-Refe rences — 1: This work: 2: iMolaro et al.l d2000h: 3: ICooke et all ll201lF) : 4: IPetitiean. Ledoux. k Sriariandl ll200Sf): 5: lEUisori et all 
pOldt); 6: iDessauees-Zavadskv et alj (120011); 7 : iPenprase et al.l ||2010|); 8: iPettini et al. I J2008I ') ■ 9: ISrianand et aL I feOlOh ; 10: 
lO'Meara et alj l|2006h : 11: IProchaska et al.1 l|2002l 'l: 12: IDessauees-Zavadskv et al] l|2003h . 
b An error estimate for this measurement was not provided by the authors. 
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Table B2. [Fe/H] and [O/Fe] in metal-poor stars 



Star Name 


T eff a 
(K) 


logg b 
(cgs) 


[Fe/H] 


[O/Fe] 


Ref. c 


riu z i yo 




4950 


1.50 


o n t 1 n i n 

—2.37 ± 0.10 


i r\ o O 1 A i n 

+0.23 ± 0.13 


1 


rliJ oOD i 




6000 


4.07 


—1.11 ± 0.06 


i / \ nn 1 1 1 n /' 

±0.29 ± 0.06 


2 


rliJ 4tv5UD 




4990 


3.04 


— 2.2i ± 0.10 


i n a *~7 1 n "i "i 

±0.47 ± 0.11 


3 


rliJ ZDlDy 




4972 


2.49 


—2.19 ± 0.10 


±0.34 ± 0.09 


3 


rliJ z i yzo 




5044 


2.67 


n n p" 1 n i n 

—2.05 ± 0.10 


i r\ o O 1 nil 

+0.22 ± 0.11 


3 


un /i coco 




5352 


3.15 


— 1.46 ± 0.10 


+0.33 ± 0.07 


3 


xiu y ( ozu 




5976 


4.16 


1 i f 1 n n 

— 1.16 ± 0.06 


i n o n i n 1 o 

+0.20 ± 0.12 


2 


IJn 1 f|QQ1 *7 




5300 


2.76 


—2.16 ± 0.10 


+0.49 ± 0.12 


3 


TJT^l 111 QQfl 

xiu iiiyou 




5694 


3.99 


— 1.03 ± 0.06 


+0.17 ± 0.07 


2 


rilJ IZZOOo 




4600 


1.10 


—2.71 ± 0.10 


+0.31 ± 0.13 


1 


rilJ IZDOo I 




4712 


1.66 


—2.76 ± 0.10 


±0.13 ± 0.14 


3 


rliJ IZODol 




5524 


4.48 


— 1.12 ± 0.06 


i n n o 1 n n o 

+0.33 ± 0.08 


2 


HD 128279 




5336 


2.95 


o i n 1 n i n 

—2.10 ± 0.10 


+0.17 ± 0.23 


3 


HD 132475 




5818 


3.95 


— 1.39 ± 0.06 


i n o o 1 n r\A 

+0.33 ± 0.09 


2 


HD 140283 




5690 


3.69 


n on i n c\r' 

—2.32 ± 0.06 


i n f n 1 n ~t n 

+0.50 ± 0.19 


2 


HD 160617 




5931 


3.77 


— 1.72 ± 0.06 


+0.22 ± 0.15 


2 


HD 166913 




6039 


4.11 


"i P" n 1 r\ c\c 

— 1.50 ± 0.06 


1 f\ O n 1 n n A 

+0.22 ± 0.24 


2 


HD 186478 




4700 


1.30 


—2.49 ± 0.10 


+0.47 ± 0.10 


1 


HD 189558 




5613 


3.91 


— 1.07 ± 0.06 


I f\ o -i 1 n n p" 

+0.31 ± 0.05 


2 


HD 205650 




5733 


4.39 


1 i o 1 n n 

— 1.12 ± 0.06 


i n o n i r\ r\c\ 

+0.30 ± 0.09 


2 


HD 213657 




6114 


3.85 


— 1.86 ± 0.06 


+0.30 ± 0.24 


2 


HD 218857 




5015 


2.78 


"i *~7n 1 n i n 

— 1.72 ± 0.10 


+0.18 ± 0.09 


3 


HD 274939 




5090 


2.79 


— 1.43 ± 0.10 


+0.35 ± 0.05 


3 


HD 298986 




6071 


4.21 


— 1.30 ± 0.06 


+0.23 ± 0.17 


2 


BD-18°5550 


4750 


1.40 


o n a I n i n 

— 2.94 ± 0.10 


i n no I ri o A 

+0.08 ± 0.24 


1 


BD-01°2582 


5072 


2.92 


i \ t \ ' y 1 n 1 n 

— 2.03 ± 0.10 


+0.31 ± 0.11 


3 


BD +17° 3248 


5250 


1.40 


1 nn i n in 

— 1.99 ± 0.10 


in ac i nil 

+0.46 ± 0.11 


1 


BD +23° 3130 


5170 


3.00 


i 1 on 1 n t \ / • 

—2.29 ± 0.06 


i n n o 1 f\ i r* 

+0.38 ± 0.15 


2 


CS 22186 - 


035 


4900 


1.50 


n Ci O 1 A 1 A 

—2.88 ± 0.10 


i n o/ -1 1 n t % /" 

+0.26 ± 0.26 


1 


CS 22873 - 


055 


4550 


0.70 


t~% q ^ 1 n ~i n 

—2.87 ± 0.10 


i n "i r\ 1 n i o 

+0.19 ± 0.13 


1 


CS 22891 - 


209 


4700 


1.00 


—3.16 ± 0.10 


±0.41 ± 0.16 


1 


CS 22892 - 


052 


4850 


1.60 


n n 1 l a i a 

—2.91 ± 0.10 


i n i a 1 n n c 

+0.14 ± 0.26 


1 


CS 22896 - 


154 


5250 


2.70 


—2.58 ± 0.10 


+0.64 ± 0.23 


1 


CS 22948 - 


066 


5100 


1.80 


—3.01 ± 0.10 


+0.54 ± 0.23 


1 


CS 22949 - 


037 


4900 


1.50 


-3.81 ±0.10 


+1.54 + 0.13 


1 


CS 22953 - 


003 


5100 


2.30 


-2.73 ±0.10 


+0.44 + 0.23 


1 


CS 22966 - 


057 


5300 


2.20 


-2.52 ±0.10 


+0.70 + 0.20 


1 


CS 22968 - 


014 


4850 


1.70 


-3.42 ±0.10 


+0.51 + 0.27 


1 


CS 29491 - 


053 


4700 


1.30 


-2.92 ±0.10 


+0.43 + 0.17 


1 


CS 29495 - 


041 


4800 


1.50 


-2.71 ±0.10 


+0.37 + 0.13 


1 


CS 29516 - 


024 


4650 


1.20 


-2.94 ±0.10 


+0.28 + 0.20 


1 


CS 29518 - 


051 


5200 


2.60 


-2.58 ±0.10 


+0.63 + 0.23 


1 


CS 30325 - 


094 


4950 


2.00 


-3.17 ±0.10 


+0.36 + 0.33 


1 


CS 31082 - 


001 


4825 


1.50 


-2.79 ±0.10 


+0.28 + 0.15 


1 



a Stellar effective temperature. 

b Surface gravity; 

References— l:ICavrel et al.l ll2004h; 2:lNissen et alj d2002h ; 3:lGarcia Perez et al.l fc006T). 
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